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systems  to  monitor  ionospheric  conditions.  The  central  question 
of  interest  is  "How  accurately  can  the  E  and  F  region  electron 
density  profiles  (ne)  be  determined  from  satellite  optical  data 
on  a  global  scale?".  We  have  addressed  this  question  as  it  1 
pertains  to  the  auroral  E  region  and  low  to  mid-latitude  daytime: 
E  and  F  regions.  Most  of  the  reported  work  is  on  the  first  of 
these  regions.  The  approach  has  been  from  first  principles  in 
which  we  start  with  either  an  incident  auroral  electron  spectrum 
or  a  solar  £UV/X-ray  spectrum  and  proceed  to  determine  the 
associated  electron  distribution  function  followed  by  its 
associated  excitation  and  ionization  rates,  optical  intensities, 
and  electron  and  ion  densities.-^- In  our  auroral  modeling,  these 
quantities  have  been  examined  unafer  a  variety  of  conditions 
defined  in  terms  of  the  incident  electron  spectrum  and  model 
atmosphere.  We  observe  some  significant  relative  variations 
among  various  UV  intensities  as  the  electron  spectrum  is 
charged.  Features  considered  have  been  01  1356  A,  N?+  3914  A, 
and  LBH  bands  at  1273  A,  1325  A,  1354  A,  and  1384  A.  The 
relative  variations  are  caused  by  the  varying  amount  of  O2  pure 
i  absorption  from  one  feature  to  the  next  and  are  observed  to  be 
|  as  much  as  a  factor  of  5  as  the  mean  energy  of  the  incident 

spectrum  is  charged  from  1  to  10  keV.  The  results  point  to  the 

I  potential  of  using  satellite  observed  UV  intensities  to  monitor 
auroral  electron  spectral  hardness.  They  point  also  to  the 
I  potential  of  monitoring  the  0  density  since  0  emissions  such 

I  as  01  1356  A  are  sensitive  to  this  density. 

In  our  dayglow  modeling,  UV  intensities  and  ne  have  been 
calculated  as  functions  of  solar  activity  and  model  atmosphere. 
Both  types  of  calculated  quantities  are  shown  to  be  sensitive 
to  solar  activity  (factors  of  2  to  3  variation  from  low  to  high 
activity).  01  emissions,  as  expected,  are  also  sensitive  to 
j  changes  in  the  0  density. 

!  Several  quantities  were  calculated  for  comparison  with 

I  data  from  the  AFGL  auroral  E  rocket  program.  These  include 
!  electron  spectra,  ion  and  electron  densities,  and  altitude  pro- 
:  files  of  several  UV  intensities.  Overall  good  agreement  is 
achieved  between  the  calculations  and  preliminary  versions  of 
data  so  far  made  available  to  us. 
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Introduction  and  Technical  Summary 


In  this  report,  we  discuss  work  undertaken  to  better  understand  the 
properties  of  optical  emissions  in  the  auroral  and  low  to  mid'latitude  daytime  iono¬ 
spheres.  The  approach  is  to  use  first-principles  techniques  to  calculate  the 
needed  quantities  under  a  variety  of  geophysical  conditions.  In  particular, 
we  have  calculated  optical  intensities  and  electron  densities  tne(z)]  noting  their 
variability  from  one  set  of  conditions  to  another.  A  long  term  goal  of  this 
work  is  to  determine  the  usefulness  of  satellite  observed  optical  intensities 
for  deducing  the  electron  density  profile.  This  goal  is  of  interest  to  as  much 
of  the  global  ionosphere  as  the  technique  proves  applicable  although  here  the 
discussion  is  limited  to  those  regions  noted  above. 

We  know  of  no  prominent  emission  features  in  the  auroral  and  daytime  iono¬ 
spheres  which  qive  direct  signatures  of  the  ions  present  unlike,  say,  the  0+  recombina¬ 
tion  emissions  in  the  tropical  nighttime  airglow  (see,  e.g.,  Chandra  et  al.,  1975  and 
Tinsley  and  Bittencourt,  1975).  0+  834  A  comes  closest  to  fulfilling  the  require¬ 

ment  since  0+  can  affect  the  834  A  intensity  through  multiple  scattering  (see 
Feldman  et  al . ,  1981  and  Kumar  et  al.,  1982).  This,  in  itself,  poses  a  difficult 
problem  which  we  shall  not  address  here.  It  then  becomes  a  two  step  process  to 
obtain  n0  from  optical  measurements.  We  must  first  obtain  a  representation  of 
the  source  (solar  ionizing  spectrum  or  incident  auroral  electron  spectrum)  fol¬ 
lowed  by  its  use  in  calculating  n  .  The  emphasis  in  this  work  is  on  how  well  the 
source  spectrum  can  be  determined  from  optical  data.  To  address  the  second  step, 
we  then  calculate  n0  examining  its  sensitivity  to  variations  in  the  source  spectrum 
and  the  major  neutral  densities. 

The  variability  of  n0  can  pose  some  serious  problems  to  the  above 
technique.  This  variability  is  caused  by  1)  time  and  spatial  changes  in  the 
source,  2)  the  same  type  of  changes  in  the  major  neutral  density  profiles,  and 
3)  mechanisms  which  produce  bulk  transport  of  the  plasma.  As  noted  above,  item  1 
shall  receive  most  of  our  attention.  Item  2  is  important  because  uncertainties 
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introduced  can  affect  both  the  accuracy  of  the  deduced  source  spectrum  and  the 
subsequently  calculated  ng  profile.  The  accuracy  of  ng  is  dependent  not  only 
on  the  assumed  source  spectrum,  but  also  on  assumed  neutral  densities  through 
chemistry  and  diffusion.  More  shall  be  said  about  item  2  in  the  next  few  oa's- 
graphs.  Item  3  has  to  due  with  winds  and  fields  and  will  be  briefly  discussed 
in  the  next  section.  Beyond  this,  however,  we  shall  restrict  ourselves  to  con¬ 
ditions  for  which  n&  is  dependent  only  on  the  ionizing  source,  chemistry,  and 
diffusion . 


The  technique  being  addressed  is  most  attractive  for  situations  in 
which  the  density  profiles  of  N2>  02,  and  0  are  known.  Then,  variability  of 
one  optical  intensity  relative  to  another  with  changes  in  either  time  or  location 
can  be  directly  related  to  variability  in  source  conditions.  The  neutral  density 
profiles  cannot  be  precisely  known,  however,  because  of  their  intrinsic  variability 
and  the  difficulty  of  precisely  describing  that  variability.  The  latter  problem 
is  due  to  the  approximate  nature  of  existing  thermospheric  density  models  (for 
recent  models, see  Jacchia,  1977;  Hedin  et  a  1.,  1977  a,b;  Barlier  et  al . ,  1978; 

Hedi r,  et  al.,  1981)  and  to  the  limitations  of  measurement  techniques.  Of  par¬ 
ticular  interest  here  are  satellite  based  measurements.  In  situ  measurements  can 
be  expected  to  be  the  most  reliable  but  may  be  insufficient  for  specifying  alti¬ 
tude  profiles  to  the  degree  of  accuracy  desired.  Optical  remote  sensing  techniques, 
on  the  ocher  hand,  do  hold  a  potential  for  measuring  column  densities  which  can  be 
used  to  scale  corresponding  volume  density  profiles.  The  next  paragraph  addresses 
this  problem.  Returning  to  the  original  issue  of  this  paragraph,  one  of  the  ob¬ 
jectives  of  the  present  work  must  be  to  determine  just  how  precisely  the  major 
neutral  densities  need  to  be  known  to  make  the  technique  under  investigation 
practical.  The  required  precision  is  dependent  on  source  conditions  and,  of 
course,  on  instrumental  conditions.  We  have  begun  to  answer  the  question  of  pre¬ 
cision  as  it  relates  to  source  conditions  which  will  be  taken  up  later  in  this 
report. 


The  problem  raised  in  the  previous  paragraph  is  whether  it  is 
possible  to  separate  the  effects  of  neutral  density  variations  from  source 
variations  when  using  optical  data.  Various  researchers  are  actively  investigat¬ 
ing  the  problem  of  sensing  the  neutral  atmosphere  (Newman  et  al.,  198?  and 
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Meier  and  Anderson,  1982)-  Their  concern  has  been  with  the  analysis  of 
satellite  observed  limb  profiles  for  features  such  as  01  1356  A  and  LBH 
bands.  We  have  concentrated  our  efforts  on  nadir  viewing  observations  and  we 
believe  that  with  the  right  choice  of  features  and  instrumental  parameters, 
the  source  and  atmospheric  variabilities  can  be  decoupled.  This  is  strengthened 
by  the  fact  that  there  are  limitations  to  the  relative  variability  between 
and  0,  densities  in  the  regions  of  peak  excitation  (for  either  auroral  or  daytime 
conditions).  Furthermore,  the  models  cited  above  are  probably  more  than  adequate 
under  most  observing  conditions  for  describing  these  relative  densities.  It  then 
becomes  attractive  to  use  selected  band  emissions  for  initially  characterizing 
the  souce  spectrum.  For  auroral  conditions,  their  relative  strengths  indicate 
the  hardness  of  the  electron  source  spectrum  through  either  of  the  basic  mechanisms 
of  pure  absorption  or  quenching,  depending  on  the  features.  Given  a  representa¬ 
tion  of  the  spectral  hardness,  the  magnitudes  of  the  observed  intensities  then 
determine  the  energy  content  of  the  spectrum.  For  daytime  conditions,  there  is 
little  variation  in  the  hardness  of  the  relevant  part  of  the  solar  ionizaing 
spectrum  and  so  here,  it  is  the  energy  content  of  the  spectrum  we  seek.  The 
magnitudes  of  the  observed  intensities  provide  this  information. 

With  an  initial  characterization  of  the  source  spectrum,  we  may  turn  to 
the  problem  of  better  determining  the  neutral  densities.  The  atomic  oxygen  density 
is  of  prime  interest  since  it  probably  experiences  the  greatest  relative  variability 
(see  above  papers  on  thermospheric  modeling)  and  has  associated  with  it  some  kpy  UV 
emission  features  (examples  are  01  1304  A,  01  1356  A  ,  and  01  2972  A).  Because 
0  is  a  minor  species  when  considered  over  the  entire  excitation  region,  its 
optical  intensities  are  sensitive  to  its  overall  density.  Thus,  given  the  source 
spectrum,  a  measure  of  say,  the  01  1356  A  intensity  (with  a  minor  contamination  from 
Ng  LBH  1354  A)  gives  a  direct  measure  of  the  magnitude  of  the  0  density  profile.  We 
may  now  repeat  the  process  to  further  improve  both  the  source  representation  and 
neutral  density  description.  Quantitative  information  related  to  the  discussion 
over  the  past  few  paragraphs  will  follow  in  subsequent  sections. 

A  summary  of  the  contract  work  now  follows.  This  will  be  given  in  terms 
of  results  and  improvements  in  models  and  codes.  Some  of  the  results  have  already 
been  presented  at  the  1981  Annual  Fall  meeting  of  the  AGU  and  in  the  papers'  by 
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Strickland  et  at.  (1982b)  submitted  to  JGR,  and  Strickland  et  al.  'I  982c;  submitted 
to  the  AIAA.  We  begin  by  listing  key  results.  Conclusions  are  based  on  tne-  assu"ptior 
that  fields  and  winds  are  not  significantly  perturbing  the  plasma.  Since  this  is 
frequently  not  the  case  in  the  F  region,  the  results  to  follow  refer  to  the  E  and 
possibly  lower  F  regions. 


1)  In  our  auroral  modeling,  selected  UV  intensities  have  been  found 

to  have  noticeable  variations  relative  to  one  another  as  the  incident 
electron  spectrum  and  model  neutral  densities  are  varied.  Some 
features  whose  intensities  show  this  behavior  are  the  LBH  bands  at 
1273  A,  1325  A,  1354  A,  and  1383  A,  Vegard-Kaplan  bands  such  as  that 
at  2672  A,  and  the  atomic  oxygen  lines  at  1356  A  and  2972  A.  The 
calculated  electron  density  profile  shows  significant  variation 
with  the  above  variations  in  the  intensities.  We  conclude  that 
good  satellite  measurements  of  intensities  such  as  the  above  should 
allow  one  to  deduce  ng  throughout  the  region  of  important  optical 
emission  to  better  than  5011.  Some  qualifications  are  called  for 
and  we  refer  the  reader  to  the  appendix. 

2)  For  daytime  conditions,  there  is  little  relative  variation  of 
intensities  as  the  solar  spectrum  changes.  Magnitudes  of  both 
intensities  and  the  np  profile  do  change,  however,  by  factors  of 
two  in  going  from  solar  minimum  to  maximum  conditions.  An  ad¬ 
vantage  of  there  being  little  relative  variation  among  intensities 

is  that  altitude  profiles  of  emissions  and  ng  remain  essentially 

constant  with  time  and  location  for  similar  solar  zenith  ancles 
and  exospheric  temperatures.  Changes  in  profile  shape  which  do 
take  place  as  a  result  of  changes  in  these  tw  parameters  cause 
no  problem  since  they  can  be  accurately  modeled.  We  conclude 
that  good  satellite  measurements  of  intensities  such  as  listed 
above  should  give  a  good  magnitude  measure  of  the  ionizing  part 
of  the  solar  spectrum  as  well  as  the  0  density  profile.  From 
this  information,  we  expect  that  in  the  E  and  F  regions  can  be 
determined  to  even  greater  accuracies  than  in  the  auroral  E  region 
under  quiescent  plasma  conditions. 
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3)  A^GL  conducted  an  auroral  E  rocket  program  in  March  of  1981  to 
observe  the  continuous  aurcra.  We  have  provided  a  variety  c*  pre¬ 
dictions  for  comparison  with  data.  Quantities  include  intensity 
altitude  profiles  of  numerous  optical  features,  ion  density  profiles, 
and  n  (z).  Some  of  the  data  have  been  made  available  to  us  in  pre¬ 
liminary  form.  We  find  good  agreement  in  the  fractional  abundances 
of  N0+,  0p+,  0+  and  N+.  Agreement  in  n&  to  better  than  50“  is  also 
obtained.  Most  of  the  calculated  intensity  profiles  alsc  compare 
favorably  in  shape  and  magnitude  with  the  data  from  the  fielded 
photometers  and  UV  spectrometer.  More  definitive  statements  will 

be  made  when  the  data  are  released  in  their  finalized  fGrns  but 
current  indications  suggest  that  the  continuous  aurca  can  be 
effectively  modeled. 

4)  A  satellite  optical  experiment  designed  to  observe  far  UV  (_■  2000  A) 
emissions  in  the  nadir  direction  will  be  limited  to  monitoring  the 
region  above  100  km.  This  is,  of  course,  desirable  for  our 
problem  of  interest  since  it  concerns  remote  sensing  of  the  ionosphere. 
The  questions  ari?ec  as  to  how  long  the  wavelength  of  it,  emission 
feature  can  be  before  Rayleigh  scattering  as  well  as  scattering  from 
clouds,  aerosols,  and  the  earth's  surface  begin  to  produce 

a  backscattered  intensity  interfering  with  the  measurement  of 
the  direct  ionospheric  component.  The  answer  depends  on  whether 
solar  radiation  is  present  or  not.  For  auroral  observations  in 
the  absence  of  sunlight,  backsca ttering  does  not  pose  a  problem 
shortward  of  3000  A.  For  daytime  conditions,  this  wavelength  is 
in  the  Schumann-Runge  band  absorption  region  near  1900  A.  A  detailed 
discussion  of  this  subject  is  given  in  Section  6. 


To  complete  the  summary,  we  note  improvements  that  have  been  made  to 
our  models  and  corresponding  codes.  As  is  usually  the  case  for  this  type  of 
activity,  the  implementations  constituted  an  important  part  of  our  overall  effort. 
The  following  list  highlights  the  most  important  of  the  improvements. 
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1)  'he  chenistry  codes  now  allow  for  diffusion  of  0^  in  the  F  region. 
Section  5  details  the  applied  algorithm. 

2)  The  electron  temperature  is  now  calculated  along  with  ion  and 
neutral  species  in  the  chemistry  codes.  Heat  conductor  .  _  ,.j, 
accounted  for  but  will  soon  be  added. 

3)  Intensities  are  now  calculated  for  many  more  features  than  ea>-lier 

in  this  program.  Most  additions  are  to  band  systems.  (Np+  IN,  N  jp, 

N?  VK,  LBH ,  and  BH). 

4)  Reaction  rate  coefficients  have  been  updated  as  data  became  available. 

1  3 

The  chemistry  of  0(  S )  and  Np  (A  1)  was  examined  in  some  detail  leadin 
to  better  descriptions  of  the  emissions  at  01  5577  A  and  01  2972  A 
and  in  the  Vegard-Kaplan  system. 

5)  Code  PEGFAC  became  functional  on  the  AFGL  system.  This  code  calculate.', 
photoelectron  spectra  and  associated  excitation  rates.  Modifications 
were  also  made  to  the  auroral  chemistry  code  for  daytime  applications. 

6)  Significant  reductions  were  made  in  disk  storage  and  runtime 
requirements  for  key  codes.  Among  the  advantages  of  such  re¬ 
ductions  has  been  much  better  “turn  around"  time  for  the  auroral 
electron  transport  code. 
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Ionospheric  Variability  due  to  Neutral  Winds 
and  Electric  Fields 


In  the  F  region  ionosphere  (above  180  km)  both  the  ion-neutral  and 
electron-neutral  collision  frequencies  are  much  less  than  the  respective  gyro- 

frequencies.  This  means  that  the  mobility  of  a  charged  particle  is  much  greater 
along  a  field  line  than  across  field  lines.  The  neutral  winds  of  the  thermo¬ 
spheric  circulation  produce  a  frictional  force  on  the  ionospheric  plasma,  but 
only  the  component  of  that  force  which  acts  along  a  field  line  is  effective  at 
imparting  motion  to  the  plasma.  Magnetic  field  lines  lie  approximately  in 
meridional  planes  and  rise  in  altitude  toward  the  equator.  During  daytime  hours, 
the  normal  thermospheric  circulation  is  from  equator  to  pole  (Bauer,  1973  and 
Evans,  1976)  which  tends  to  drive  the  plasma  down.  This  reduces  not  only  the 
altitude  of  the  ionization  peak,  but  also  the  peak  concentration  since  chemical 
losses  are  larger  at  the  lower  altitudes.  At  middle  latitudes  the  effects  of 
these  winds  can  be  as  important  as  diffusion,  resulting  in  a  lowering  of  the 
\ 2  peak  by  ^  50  km  and  a  reduction  in  peak  electron  density  by  n-  35c.-  (Hargreaves, 

1979  and  Banks  and  Kockarts,  1973B).  However,  the  thermospheric  circulation  is 
quite  variable  on  a  daily  basis  which  makes  the  ionospheric  variability  difficult 
to  predict.  One  source  of  variability  is  the  energy  deposited  at  high  latitudes 
during  geomagnetic  substorms  which  can  alter  the  global  thermospheric  ciiCuldci. 

0 roV0rsinn  th«  direction  of  flow  (Evans,  1976  and  Prolss,  1983).  Fuithermore, 
the  resistance  of  ions  to  cross  field  motion  can  alter  the  thermospheric  winds  in  the 
altitude  r^nime  above  300  km,  ^ur.thnr  comnlicatinn  the  nicturp. 

In  the  E  region  ionosphere  (105-160  km)  the  ion  cross  field  mobility 
is  considerably  larger  than  that  of  electrons.  In  this  altitude  regime  neutral 
winds  induce  a  charge  separation  which  results  in  a  global  electric  field  system. 
Because  magnetic  field  lines  are  nearly  equipotential s,  this  electric  field 
system  is  mapped  into  the  F  region  where  it  causes  plasma  (E  x  B)  drifts (Har¬ 
graves,  1979).  One  may  think  of  the  E  region  winds  acting  as  a  dynamo  genera¬ 
tor  while  the  F  region  plasma  responds  as  an  electric  motor.  The  effect  of 
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the  electric  field  is  generally  smaller  than  that  due  to  F  region  winds 
and  varies  throughout  the  day. 

In  addition  to  the  large  scale  effects  caused  by  thermospheri c  winds 
and  the  dynamo  electric  field,  there  are  small  scale  structures  -  ionospheric 
irregularities  -  which  have  a  variety  of  sources.  The  irregularities  are 
present  in  both  the  £  region  and  F  region  and  at  all  latitudes  (Ossakow,  1979 
and  Fejer  and  Kelley,  1980).  At  auroral  and  equatorial  latitudes  the  primary 
causes  of  ionospheric  irregularities  are  plasma  instabilities,  although  the 
relative  importance  of  neutral  dynamics  in  the  explanation  of  equatorial  spread 
F  is  subject  to  debate.  At  mid-latitudes  the  situation  is  less  well  understood. 
Wind  shears  are  known  to  produce  sporadic  E  layers  while  traveling  ionospheric 
disturbances  have  been  explained  by  gravity  waves  propagating  from  below  or 
from  the  auroral  zone  (Hargreaves,  1979).  A  number  of  plasma  mechanisms 
have  been  proposed  to  explain  smaller  scale  structures,  but  verification  of 
tnese  hypotheses  awaits  more  complete  sets  of  observational  data  (Fejer  and 
Kelley,  1980). 

Both  E  region  and  F  region  thermospheric  winds  can  have  a  significant 
influence  on  the  electron  distribution  above  200  km  or  so.  They  are  hard  to 
observe  and  continue  to  be  the  subject  of  theoretical  and  experimental  research. 

If  the  winds  are  known,  their  effects  on  the  ionosphere  can  be  determined  (Banks 
and  Kockarts,  1973B).  The  principal  implication  for  satellite  monitoring  of  the 
ionosphere  is  that  some  method  of  measuring  or  calculating  thermospheric  winds 
must  be  developed.  Global  models  of  the  thermospheric  circulation  provide  a  usefu1 
beginning  but  the  effects  of  geomagnetic  activity  must  also  be  included  if  reason¬ 
able  accuracy  is  to  be  obtained.  Since  ionospheric  irregularities  are  less  well 
understood,  particularly  at  mid-latitudes,  they  present  a  formidable  problem 
for  satellite  monitoring.  These  problems  are  not  addressed  further  in  this  report, 
but  remain  important  research  subjects. 
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Section  v 

Computer  Code  Development 

Substantial  improvements  were  made  in  key  computer  codes  during  the 
course  of  this  work.  These  improvements  will  be  noted  below  as  we  briefly  de¬ 
scribe  the  function  of  each  code  which  has  been  applied  to  our  auroral  and  day¬ 
time  modeling  efforts.  The  following  list  identifies  these  codes: 

1)  MX 

2)  B3C 

3)  PRATES 

4)  CHEM1 

5)  PEGFAC 

and  6)  CHEM2 

The  first  four  are  used  to  model  electron  transport,  chemistry,  and 
photon  emission  in  the  auroral  ionosphere  and  come  from  previous  work,  part  of 
which  has  been  recently  sponsored  by  AFGL  (Strickland,  1981).  The  remaining 
two  codes  are  used  for  daytime  ionospheric  modeling.  PEGFAC  conies  from  indepen¬ 
dent  work  while  CHEM2  comes  from  modifications  of  CHEM1  performed  under  this  con¬ 
tract.  Information  on  individual  codes  now  follows. 

Code  MX  generates  matrices  approximating  the  Boltzmann  collision 
integral  for  electron  scattering  and  energy  loss  as  well  as  secondary  electron 
production.  A  description  of  the  applied  transport  equation  may  be  found  in  the 
appendix.  Two  improvements  to  MX  have  recently  been  made.  One  of  these  decreases 
the  running  time  by  more  than  a  factor  of  10.  This  was  made  possible  by  replacing 
a  numerically  based  matrix  inversion  technique  by  an  analytic  one.  The  matrices  to 
be  inverted  are  of  dimension  6x6  and  can  be  decomposed  into  a  product  of  matrices 
with  dimensions  3x3  and  2x2.  This  enables  one  to  apply  analytic  techniques  to 
the  inversion  without  being  overly  burdened  with  algebra.  The  second  improvement 
has  led  to  a  significant  reduction  in  the  number  of  matrix  elements  to  be  stored 
on  disk.  Basically,  redundent  elements  have  been  removed  before  writing  to  disk 
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and  re-inserted  at  the  time  they  are  needed  to  perform  the  transport  calculation. 


Code  B3C  is  the  Boltzmann  three  (3)  constituent  electron  transport 
code,  (see  Strickland  et  al.,  1976  and  the  appendix).  It  calculates  the  pi*- 
ferential  electron  spectrum  o(z  ,£,p)  in  units  of  electrons/cm  -s-eV-sr  where 
z,  E,  and  p  refer  to  altitude,  energy,  and  cosine  of  the  pitch  angle.  B3C 
currently  allows  for  35  altitudes,  50  energies,  and  20  p  values.  The  most 
noteworthy  improvement  recently  made  allows  for  local  energy  deposition  at 
low  energies,  thereby  significantly  reducing  memory  requirements  and  running 
time.  The  local  energy  deposition  approximation  has  been  found  to  be  valid  so 
long  as  most  of  the  energy  content  of  the  incident  spectrum  lies  above  the 
transition  energy  chosen  for  going  from  the  transport  to  the  local  description. 
For  an  incident  spectrum  characterized  by*  say,  a  2  keV  Maxwellian  distribution, 
the  transition  energy  can  be  chosen  to  be  as  high  as  several  hundred  eV. 

Code  PRATES  calculates  a  variety  of  volume  excitation  and  ionization 
rates  using  <p(z  ,E,p).  The  rates  provide  the  driver  for  chemistry  modeling  in 
CHEM1 .  PRATES  is  now  significantly  larger  than  it  was  earlier  in  this  work 
following  the  above  described  modifications  to  B3C.  The  expansion  is  due  to 
the  transferal  of  numerous  routines  from  B3C  which  performed  various  operations 
following  the  actual  transport  calculations  giving  p(z  , E,u).  One  such  operation 
is  the  testing  of  energy  conservation.  The  reason  for  the  above  transferal  was 
to  increase  the  "turn  around  time"  of  B3C  runs  which  had  become  approximately 
one  day. 

Code  CHEM1  solves  a  set  of  time  dependent  rate  equations  for  numerous  ion  and 
neutral  species.  It  currently  also  calculates  optical  intensities  for  a  large 
number  of  atomic  and  molecular  emission  features.  CHEM1  has  been  undergoing 
continued  development  throughout  the  given  contract  period  of  performance. 

Some  of  this  effort  has  involved  updating  chemical  reaction  rate  coefficients 
and  adding  emission  parameters  for  additional  optical  features.  The  rest  has 
led  to  CHEMl's  ability  to  treat  diffusion  and  calculate  the  electron  temperature, 
Tg.  A  discussion  of  the  diffusion  work  appears  in  Section  5.  Tg  is  currently 
calculated  assuming  no  conduction  but  our  plans  call  for  the  removal  of  this 
restriction  during  the  next  phase  of  work  on  the  code. 
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Code  PEGFAC  calculates  the  differential  photoelectron  spectrum 
0(z  ,E)  in  units  of  electrons/cm  -s-eV-4u  sr.  Section  8  details  PEGFAC's 
history,  computational  model,  input  data,  and  selected  results.  Most  of 
the  required  effort  related  to  PEGFAC  was  directed  to  its  transferal  to  the 
AFGL  CDC  computing  system  from  a  non-CDC  system  upon  which  it  had  resided. 

Code  CHEM2  serves  in  the  same  capacity  for  daytime  ionospheric 
modeling  as  CHEM1  does  for  auroral  modeling.  As  noted  above,  CHEM2  was  de¬ 
veloped  using  CHEM1  as  its  basis.  Differences  between  these  codes  are  minor 
as  of  this  writing  and  primarily  relate  to  the  handling  of  input  data  generated 
by  respective  electron  flux/excitation  rate  codes. 


't'i  c  n 


The  Chemistry  of  Q(*S)  and  (A  3Z) 


Because  the  0(^S)  emission  feature  at  5577  A  is  one  of  the  most 
prominant  auroral  features,  it  should  be  included  in  any  comprehensive  auroral 
model.  Although  non-ionospheric  emission,  scattering,  and  albedo  effects  pre¬ 
vent  the  useful  nadir  observation  of  this  emission  from  satellites,  its  calcula¬ 
tion  does  provide  a  useful  check  on  the  chemistry  model.  In  addition,  the 
chemistry  of  0(*S)  is  coupled  to  that  of  Ng(A  3Iu+)  which  gives  rise  to  the 
Vegard-Kaplan  bands.  These  bands  may  prove  useful  for  monitoring  the  hardness 
of  precipitating  electron  spectra  from  satellites.  Consequently,  we  have  di¬ 
rected  some  of  our  effort  to  include  the  best  available  information  on  0(^S) 
and  Np(A)  in  our  model . 


4.1  The  N^A)  State  and  the  Vegard-Kaplan  System 

The  N^(A  3Iy  )  metastable  state  is  populated  by  direct  excitation 

tr-om  the  ground  state  and  by  cascade  from  higher  lying  triplet  states.  Since 

the  cascade  occurs  in  a  much  shorter  time  than  the  lifetime  of  the  A3  I  state, 

we  can  calculate  the  total  production  rate  for  N2(A  3Zu+)  by  summing  the  electron 

impact  cross  sections  of  all  the  trip’et  states.  Since  the  cross  sections  for 

tne  B1  ::u  ,  L  rg+,  and  D  cross  sections  are  much  smaller  than  the  others, 

we  have  neglected  them  in  the  sum.  Also,  because  approximately  half  of  the  C  1 

'  u 

state  dissociates,  we  only  include  half  of  its  cross  section: 


’total  ■  i<A  V>  *  o(B  ll  g)  *  ]  o(c  \)  *  „(W  \) 


Following  Sharp  and  Torr  (1979)  we  concentrate  on  the  v1  =  1  vibrational 
level.  Using  the  relative  populations  of  the  different  vibrational  levels  as 
calculated  by  Cartwright  (1973),  and  the  transition  probabi 1 i ties  of  Shemansky 
(1969)  we  have  calculated  the  emission  rates  for  each  vibrational  level  assuming 
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no  quenching.  In  the  steady  state  and  in  the  absence  of  quenching,  the  pro¬ 
duction  rate  must  equal  the  emission  rate.  Using  our  calculated  emission  rates, 
we  have  calculated  the  relative  production  rates  of  the  v'  -  0  and  v'  -  1  levels, 
and  we  assume  that  these  rates  are  valid  even  when  quenching  becomes  important. 


There  are  only  two  loss  mechanisms  for  the  N2(A  .  )  state:  quenching 

and  emission.  Except  at  the  lowest  altitudes  where  Oj  becomes  important,  the 
dominant  quencher  is  atomic  oxygen.  For  completeness,  we  have  also  included 
quenching  by  and  NO.  The  density  of  any  vibrational  level  is  determined  by 
the  balance  between  production  and  loss: 


P(A,v' )  =  (qy,  +  Ay, v„)  n  (A ,v'  ) 


where  qy,  =  l  kj(v')n^  is  the  sum  of  the  quenching  rates  due  to  each  neutral 
species,  k^(v’)is  the  quenching  coefficient  for  the  v'  level,  and  Av,v„ 
is  the  transition  probability  from  the  v'  level  of  the  A  state  to  the  v"  level 
of  the  ground  state.  The  actual  emission  from  any  vibrational  level  is 


e(A,v'  )  =  (I  A  ,  ,,)  n  (A,v’ ) 
v "  v  v 


Eliminating  the  concentration  of  the  v'  level  between  Equation  (1)  and  (2) 
results  in 


Av.p(A,v') 

■rrnrr 

v 
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where 


We  have  adopted  the  transition  probabil ities  given  by  Vallance 
Jones  (1974).  Since  these  are  somewhat  larger  than  those  reported  by  ur:.;.'  gnr 
(1978),  they  may  need  to  be  revised  in  future  calculations.  The  quenching  co¬ 
efficients,  k ^ ( v ' )  were  taken  from  the  review  by  Torr  and  Torr  (1982).  Unfortunately, 
the  quenching  coefficient  for  the  most  important  species,  0,  is  also  the  most 
uncertain.  A  recent  laboratory  measurement  (Piper  et  al .  ,  1981)  found 
kQ(v'  =  0)  =  2.8  x  10’11  and  kg ( v '  =  1)  =  3.4  x  10’11  cm^s"1.  This  is  in  general 
agreement  with  the  values  adopted  by  Cartwright  (1978)  based  on  the  ground  based 
observations  of  Vallance  Jones  and  Gattinger  (1976).  However,  these  values  are 
in  sharp  disagreement  with  the  rocket  determination  by  Sharp  (1971)  who  found 
Kq ( v  =  0)  =  9  x  10’11  and  k0(v'  =  1)  =  2  x  10"10.  Until  this  discrepancy  is 
resolved,  any  calculation  of  the  Vegard'Kaplan  emissions  is  subject  to  some  un¬ 
certainty. 


Using  the  laboratory  values  of  kg(v'),  we  have  calculated  the  Vegard- 
Kaplan  emission  for  several  incident  electron  energies.  Volume  production,  quench¬ 
ing,  and  emission  rates  for  1  keV  and  5  keV  Maxwellian  cases  are  displayed  in 
Figures  1  and  2,  respectively.  Column  emission  rates  are  shown  in  Figure  3  as  a 
function  of  incident  electron  energy.  The  total  emission  from  each  vibrational 
level  and  the  two  most  prominant  bands  of  each  progression  are  displayed.  Also 
shewn  is  the  energy  dependence  of  the  total  emission  when  the  larger  quenching 
coefficients  of  Sharp  (1971)  are  used.  Some  noticeable  differences  in  emission 
are  seen  to  occur  when  using  the  Sharp  coefficient  values  in  place  of  the  Piper 
values.  Since  the  weight  of  evidence  favors  the  Piper  et  al.  (1981)  quenching 
coefficients,  it  is  likely  that  we  can  expect  a  factor  of  3  to  5  variation  in 
some  of  the  stronger  bands  as  the  incident  electron  energy  changes  from  1  to  5  keV. 
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Altitude  (km) 


PRODUCTION  AND  LOSS  RATES  (cm'V1) 


Figure  1  Volume  Production,  Emission,  and  Loss  Rates  for  the  v'  =  0 

Progression  of  the  Vegard-Kaplan  System.  The  Incident  Electron 
Spectrum  was  a  1  keV  Maxwellian.  We  have  used  the  quenching 
coefficients  of  Piper  et  al.  (1981)  and  the  Jacchia  (1977)  model 
atmosphere. 
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Altitude  (km) 


MAXWELLIAN  CHARACTERISTIC  ENERGY  EQ(keV) 


Figure  3  Energy  Dependence  of  Veqard-Kaplan  Emission  for  an  Energy 

-2-1 

Deposition  Rate  of  1.0  erg  cm  sec  in  the  Jacchia  (1977) 
Model  Atmosphere.  Total  Emission  From  the  v'  =  0  and  v'  =  1 
Progressions  are  Shown  for  two  Different  Values  of  the  Atomic 
Oxygen  Quenching  Coefficient  (see  text).  In  Addition^  the 
Emissions  Rates  of  the  two  Most  Prominent  Bands  in  Each  Pro- 
aression  are  Shown  for  the  Most  Probable  Value  of  the  Quench 
ina  Coefficient. 


4.2  The  Chemistry  of  0(^S) 

Aside  from  direct  excitation  by  energetic-electrons,  the  principa1 
sources  of  0(*S)  are  (Torr  and  Torr,  1982): 


02+  +  e  -  0  +  0(1S) 

n2(a3eu+)  +  0  -  n2  +  0(1S) 

02+  +  N  -  N0+  +  0  (!S) 


The  temperature  dependent  reaction  rate  for  dissociative  recombination  of 
02+  has  been  given  by  Roble  and  Rees  (1977)  which  is  similar  to  the  rate  de¬ 
termined  from  atmospheric  explorer  measurements  (Torr  and  Torr,  1978).  The 
branching  ratio  varies  from  2%  to  102  depending  on  the  vibrational  level  of  02+ 

(Torr  and  Torr,  1982).  We  have  adopted  a  value  of  8%.  For  the  third  reaction,  0?+ 

-113-1  c 

with  N,  we  have  adopted  a  value  of  2.5  x  10  cm  s  which  falls  in  the  middle  of 

of  the  range  reported  by  Torr  and  Torr  (1982). 

Sharp  and  Torr  (1979)  determined  the  rate  of  production  of  0(^S; 
from  the  quenching  of  N2(A  )  by  0.  Using  the  population  of  the  v'  =  1  level 

of  N,,(A  +  ),  whose  calculation  is  described  in  the  preceeding  section,  and  an 

c  U  -11  3  -1 

efficiency  of  25%,  they  obtained  the  effective  rate  3.6  x  10  cm  s  .  Note  that 

this  rate,  when  multiplied  by  the  atomic  oxygen  concentration  and  the  population 

of  the  v1  •-  1  level  of  N2(A  Iu  ),  gives  the  total  production  of  0(  S)  from  a1 1 

vibrational  levels  of  N2(A). 

The  principal  chemical  losses  we  have  included  are  quenching  by  0,  NO, 
and  02-  Of  these,  quenching  by  02  is  the  most  important  one.  We  have  used  a 
temperature  dependent  coefficient  for  this  reaction  from  Hyman  and  Julienne  (1975). 
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We  are  now  calculating  the  0(^S)  density  and  its  corresponding 
5677  A  emission  using  the  above  described  information.  Figure  4  gives  an 
example  of  a  zenith  viewing  5577  A  column  emission  rate  for  a  incident  elec- 
tron  spectrum  given  by  a  1  keV  Maxwellian  with  an  energy  content  of  1  erg/cm  - 
Also  included  is  the  emission  for  N^+  at  3914  A  since  it  is  commonly  measured 
along  with  5577  A  in  the  same  experiment.  The  observed  ratio  for  the 
column  emissions  (5577  A/ 3914  A)  is  generally  between  1  and  2  which  provides 
us  with  some  confidence  that  we  are  achieving  a  reasonable  description  of 
the  5577  A  emission. 


Altitude  (km) 


4tI  (kR) 


Figure  A  Zenith  Viewing  Column  Emission  Rates  for  01  5577  A  and 

Ng  3914  A.  The  Incident  Electron  Spectrum  is  Character¬ 
ized  by  a  Maxwellian  Distribution  With  Characteristic 
Energy  and  Energy  Content  as  Shown. 
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The  Diffusion  of  0+  in  the  F  Region 

In  the  dayside  and  auroral  ionospheres,  the  F2  region  is  predominantly 
0+(^S)  (hereafter  referred  to  as  0+).  Calculations  of  the  0+  density  (which  is 
equivalent  to  the  electron  density)  need  take  into  account  only  those  processes 
which  involve  that  particular  species.  This  results  in  considerable  simplifica¬ 
tion  of  the  diffusion  eauation,  which  may  be  written  as. (Banks  and  Kockarts, 

1973  B.  d.  170)  as  : 

=  q  -  In  -  rj  (d)  (3a 

:  =  nv  =  -  Da  sin2  I  +(j-  jf  +  p-  )  n  (3b 

with  the  following  definitions  : 

-f  O 

i  =  0  flux  in  ions/cm  -s 

n  =  0+  concentration 

v  =  z  component  of  ion  (and  electron)  velocity 

a  =  production  rate 

L  =  loss  rate  per  0+  ion  due  to  chemistry  (loss  by  recombination, 

2 

given  by  an  ,  is  negligible) 

Tp  =  plasma  temperature  (ion  temperature  +  electron  temperature) 

Hp  =  plasma  scale  height  (kT  /m^g) 

Da  =  ambipolar  diffusion  coefficient  (kTp/m.vin) 
vin  =  ion  neutral  collision  frequency 

I  =  magnetic  dip  angle,  measured  from  horizontal  and  positive  in 
the  northern  hemisphere 

?! 


in  writing  Equation  (3)  we  have  neglected  inertial  terms,  cross-fit'o  ioi 
motion  ,  and  neutral  winds. 


When  the  expression  for  the  flux  (3b)  is  substituted  into  the  con¬ 
tinuity  Equation  (3a) ,  the  result  is  a  linear  parabolic  partial  differential 
equation : 


2  o  ,  3D  -i  ,  3T 

on  _  n  •  T  o  n  _  •  _  £  t  r  1  ^  *  •  1  P  t  c?  n 

-t  -  Da  sin  I  — 7  +  Da  sin  *-0  3z“  +  H~  +  T  3z  ^  37 
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+  r 
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-nr  )3 


t  :n 


(4) 


In  general  the  coefficients  can  depend  on  both  altitude  and  time  so  that 
numerical  methods  must  be  used  to  obtain  a  solution.  We  are  interested  in 
the  steady  state  solution  which  results  from  time  independent  coefficients. 
iNevertne  1  ess ,  the  simplest  method  for  obtaining  the  steady  state  solution  is 
to  choose  a  reasonable  initial  condition  and  then  integrate  Equation  (4)  un¬ 
til  a  steady  state  is  reached.  This  has  the  additional  advantage  that  the 
same  computer  code  may  be  used  for  problems  with  time  varying  coefficients 
and  boundary  conditions. 

To  obtain  the  solution  we  must  specify  not  only  an  initial  condition, 
but  also  boundary  conditions.  The  lower  boundary  is  specified  at  z  and  is  made 
low  enough  that  chemical  equilibrium  prevails.  The  resulting  boundary  con¬ 
dition  is 


n(z0)  =  q(zo)/L(zo) 


(51 


The  altitude,  zn  refers  to  the  upper  boundary  and  must  be  chosen  to  be  above  the 
electron  density  peak  but  below  the  altitude  at  which  H+  becomes  the  dominant 
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ion.  Following  Strobel  and  McEl roy  (1970)  we  have  choosen  the  upper  bound arv 
condition  to  be 


q(z)dz  n,  -  q(zn)H0(zn) 


(6) 


where  Hq(z  )  is  the  scale  height  of  neutral  oxygen  at  z^.  Using  Equation  (3b), 
the  upper  boundary  condition  takes  the  form 


F(z  ) 
n  oZ 


+  G(z  )n(z  )  =  /  q(z)dz 
n  n  / 


with  F(z)  =  Da  sin2  I  and  G  -  Dg  sin2  J  [j-  +  ^-)  , 

P  P 

To  obtain  an  approximate  numerical  solution  of  Equation  (4)  we  have 
adopted  the  Crank-Ni col  son  scheme  proposed  by  Ames  (1977  ,  o.  67).  First  we  re¬ 
write  Equation  (4)  as 


~  =  a ( 2 )  It  +  b(z)  +  c{z)n  +  q(z,t) 

9z 


where  a[z)  =  D  sin  I 

a 


,  3D  ,  ,  3T 

b(z)  =  a(z)  [p-  -gj-  +  fl-  +  r-  -gf  3 

a  p  p 


■(zl  -  a ( z  1  T  1  (  1  +  1  9TP )  +  2_  (i—  +  1—  — )]  -  L 

"  z )  "  [TT  IT  +  T~  Tz  3z  Hp  Tp  3z 

a  n  n  r  r 


p  p 


23 


and  we  have  assumed  that  a,b,  and  c  are  independent  of  time.  We  choose 
an  altitude  grid  of  n  +  1  points  numbered  from  zero  at  the  lower  boundary, 
z  ,  to  n  at  the  upper  boundary,  zn>  The  grid  spacing  is  uniform  and  given  by 
(zn  -  z  )/n.  The  time  step  is  chosen  to  be  At  and  the  following  notation 
is  adopted: 


2.  =  z  +  i:.z 
1  0 


*  j-lt 


"u  =n(2i'  V 


Vij  =  ni+l/2,j  *  n i - 1  / 2 , j 


ux"ij  =  1  (Vl/2,j  +  ni-l/2,j) 


The  discrete  version  of  (8)  is  (for  0  <  i  <  n) 


“t  (ni,j+l  -  nij}  s  ^7  (ai  {x  +  Az  Wx  +  Az2  ci)(ni,j+l  +  nij} 


+  2  (qiJ+l  +  qij} 


(10 


For  i  =  0,  the  equation  is 


o  ,j+l  '  C 


,j  +  l 

o,  j+1 


(n 


where  lQ  j  =  L(zo>tj)  is  the  chemical  loss  term.  For  i  =  n,  the  equation 
is  the  same  as  Equation  (10)  except  that  we  must  replace  quantities  evaluated  at 
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i  =  n  +  1  using  the  boundary  condi tion  (7 ) : 


=  n 


n+1 ,  j  n-l,j 


Gn 

r~ " 

n 


n,j 


(12) 


where  =  F(zn),  F„  =  G(z  )  and  Q 
n  n '  n  n 


-‘I 


q(z)dz 


When  the  indicated  operations  and  substitutions  are 
unknown  quantities  (n^^,  n.>j+1,  n.+1J+1)  are 
side.  Equations  (10  -  12)  may  be  written  in  matrix 


carried  out  and  the 
isolated  on  the  left 

form  as 


1  "j  *  Cj 


(u) 


where  n.  is  a  vector  whose  components  are  the  values  of  the  density  at  each 

J 

point  of  the  altitude  grid  at  time  t,.  The  matrices  A  and  B  are  tridiagonal 

J 

matrices  whose  non-zero  components  are 
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The  solution  of  Equation (13)  is  simplest  when  A  is  positive  definite,  or 
at  least  diagonally  dominant  and  irreducible.  Jrreducibil ity  is  assured 
if  A .  .+j  =  f •  t  0  for  i  =  1,2,3. ..n-1.  Diagonal  dominance  requires 
t  A  i ,  i  J  +  |A^  for  all  i.  This  condition  can  be  met  by 

requiring 


2a  i 

M  <  min  (r— ) 
i  i 


At  < 


In  the  auroral  ionosphere,  neither  condition  is  very  restrictive:  Az  <  2b  km 
and  At  <  40  sec  are  generally  sufficient.  Of  course,  stability  and  con¬ 
vergence  considerations  place  further  restrictions  on  A z  and  At. 

With  Az  and  At  chosen  so  that  the  matrix  A  is  diagonally  dominant. 
Equation  (13)  may  be  solved  using'Gaussian  elimination  without  pivoting. 
Furthermore,  when  the  coefficients  a,b,  and  c  are  independent  of  time, 

A  can  be  factored  once  at  the  beginnning  and  the  resulting  decompo¬ 
sition  used  at  each  time  step.  A  computer  code  incorporating  these 
features  has  been  written  and  tested  and  has  been  integrated  into  the 
existing  chemistry  code.  Preliminary  results  indicate  that  it  increases 
the  execution  ume  of  the  chemistry  code  by  only  a  small  amount.  An  example 
of  the  results  is  shown  in  Figure  5.  Although  chemical  equilibrium  is 
established  quickly  in  the  E-region,  diffusive  equilibirum  is  reached  only 
after  several  hours.  Application  of  this  computer  code  to  the  daytime 
ionosphere  is  described  in  Section  8.5. 
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Tne  Effect  of  Rayleigh  Scattering  on  Nadir  Observed  UV  Intensities 
From  Satel  1  i  tes 


6.1  General  Comments 

We  wish  to  address  the  problem  of  possible  contamination  of 
ionospherically  produced  UV  intensities  due  to  the  backscatterina  of  either 
solar  or  auroral  produced  radiation  by  Rayleigh  scattering.  Viewing  is  to  be 
downward  (not  necessarily  true  nadir)  from  satellite  altitudes.  The  wave¬ 
length  region  to  be  addressed  is  2000  to  4000  A.  Below  2000  A,  there  is  nc 
problem  caused  by  Rayleigh  scattering  due  to  the  strength  of  pure  absorption 
by  02.  Longward  of  3000  A,  Rayleigh  scattering  is  not  the  only  possible 
source  of  contamination'.  Scattering  from  clouds,  aerosols,  and  the  earth's 
surface  become  possible.  Here,  we  shall  restrict  ourselves  to  Rayleigh 
scattering  but  must  keep  in  mind  these  other  mechanisms  when  dealing  with 
features  such  as  N^+  3914  A. 

We  shall  see  that  there  can  be  serious  contamination  from  back- 
scattered  solar  radiation.  The  key  to  this  problem  is  both  the  strength 
of  the  solar  flux  and  the  width  of  the  instrumental  bandpass.  Let  us 
consider,  for  example,  a  spectrometer  with  a  resolution  of  1  A.  The 
width  of  the  emission  feature  to  be  measured  will  be  a  small  fraction  of 
l'A.  Backscattered  solar  radiation  mixed  with  this  intensity,  on  the  other 
hand,  will  cover  the  entire  bandpass  width.  This  will  place  more  severe 
restrictions  on  observable  wavelength  regions  for  daytime  conditions  in  contrast 
to  nightime  (auroral  included)  conditions  where  the  contaminated  intensity  has 
the  same  line  shape  characteristics  as  the  intensity  of  interest. 

D.E.  Anderson  Jr.  of  the  Naval  Research  Laboratory  kindly  provided 
some  of  his  time  to  help  us  begin  this  analysis.  He  also  has  provided  the 
computer  generated  plots  to  follow.  The  multiple  scattering  effects  to  be 
discussed  are  based  on  information  in  two  of  these  plots. 
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6.2 


Basic  Expressions 


The  quantity  of  interest  is  the  intensity  I  in  units  of  photons 

7  2 

(ph)/cm  -s-sr-A  for  a  wavelength  distributed  source  or  ph/cm  -s-sr  for  a 

discrete  source.  The  equation  for  I  is 


I(u)  =  Ij(u)  +  IR(-) 


(14) 


where  p  is  the  direction  cosine  relative  to  the  downward  vertical  direction 
(here,  we  are  interested  in  p  >  0),  Ij  is  the  intensity  component  due  to 
direct  emission  from  the  ionosphere,  and  IR  is  the  intensity  component  from 
below  the  ionosphere  due  to  Rayleigh  scattering.  In  the  discussion  to  follow, 
we  assume  a  plane-parallel  geometry  and  isotropic  scattering  (justified  in  this 

subsection).  Then,  either  component  has  the  following  integral  form: 

(15) 

where  t  is  the  optical  depth  for  self  absorption,  S(t)  is  the  source 

2 

function  in  units  of  ph/cm  -unit  depth  t  -  s,  and  T  is  a  transmission  func¬ 
tion  accounting  for  self  and  pure  absorption.  For  a  continuum  source,  the 
unit  A  1  must  be  added  to  the  units  of  S.  The  optical  depth  in  differential 
form,  is 


Ku)  = 


-f 


S(jriT(iA0  dt/y 
4tt 


di  =  ondz  (16) 

2 

where  o  is  the  self-absorption  cross  section  in  cm  and  n  is  the  density 
of  the  sel f-absorbers. 
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The  source  function  S( • )  is  related  to  the  total  volume  emission 

rate  by 


Si  ■  )  *-  S(z)/n 


(17 


where  the  rate  S(z)  has  units  of  ph/cmJ-s.  We  choose  to  use  S  for  either 
the  source  function  or  total  volume  emission  rate  and  let  the  argument 
distinguish  between  them.  We  introduce  S(z)  since  it  is  the  quantity  to  use  if 
one  is  dealing  with  an  optically  thin  emission.  We  will  be  considering  wave¬ 
lengths  which  apply  to  both  optically  thick  and  thin  media.  Equations  (16)  and 
(17)  allow  us  to  express  Equation  (15)  in  its  alternate  form,  namely 


K. 


/ 


S(z)  T(z,u)  dz/u 
4- 


(18 


where  we  have  followed  the  same  convention  with  T  as  with  S,  namely,  in  using 
z  and  :  as  arguments  of  the  same  function.  Equation  (18)  is  still  general,  in 
that  it  applies  to  either  optically  thick  or  thin  media.  For  the  latter  case, 
it  becomes 


l(w)  =  f  So  (z)e't/;'  dz/u  ;i° 

*'  4m 

where  5  (z)  is  the  single  scattering  or  initia1  volume  em  sston  rate  and 
t  is  the  optical  depth  for  pure  absorption  which  may  or  may  not  be  zero. 

Similar  to  S,  SQ  is  expressed  as  a  function  of  either  i  or  z. 

The  relationship  is  that  given  in  Equation  (17).  There  are  several  possible 
mechanisms  for  producing  SQ.  Those  of  interest  to  us  in  the  ionosphere  are 
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1)  auroral  electron  impact  excitation 

2)  photoelectron  impact  excitation 

3)  direct  excitation  by  solar  photons. 


and 

Those  of  interest  to  us  below  the  ionosphere  are  the  Rayleigh  scattering  of 

1)  auroral  produced  radiation 
ancl  2)  solar  radiation. 

Another  mechanism,  excitation  by  chemical  reaction s,  is  not  of  interest  here. 

If  the  medium  is  optically  thick,  then  S  will  be  greater  than  S  . 

0 

This  is  due  to  radiation  imprisonment  caused  by  the  multiple  scattering  of 
the  photons  of  interest  In  tins  case,  the  relationship  between  S  and 
may  be  expressed  through  the  integral  equation  of  radiative  transfer.  It 
i  s 


S(i)  ■  s„t- 


/ 


S(t')  H( ; 1 


)  di  ' 


(20) 


where  H  is  a  kernal  giving  the  probability  of  photons  emitted  at  ;'  reaching 
level  t  and  being  scattered  there.  Detailed  discussions  of  properties  of 
this  equation  for  resonance  line  radiations  are  given  by  Strickland  and 
Donahue  (1970)  and  Strickland  and  Rees  (1974).  Anderson  and  Meier  (1979) 
have  applied  this  equation  to  the  problem  of  Rayleigh  scattering.  They  observed 

that  adequate  solutions  to  Equation  (20)  as  well  as  Eouation  (14)  can  be  obtained 
assuming  isotropic  scattering  in  place  of  the  more  complicated  scattering 
given  by  the  Rayleigh  phase  function.  This  was  determined  by  making  com¬ 
parisons  with  'exact'  solutions  obtained  by  the  Monte  Carlo  method.  Later 

we  will  show  examples  of  S  obtained  by  Anderson  for  daytime  conditions  over 
the  wavelength  range  from  2000  to  4500  A. 
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6.3  Derivation  of  the  Intensity  Equation 


It  is  instructive  to  relate  Equation  (14)  to  what  is  actually 
observed  by  a  detector.  Consider  to  detector  to  be  situated  at  satellite 
altitudes  viewing  the  ionosphere  and  lower  atmosphere  at  some  angle 
(cos  *  *)  relative  to  the  downward  vertical  direction.  Figure  6  depicts 
the  observing  conditions.  Let  5N  be  the  rate  of  photons  received  from 
the  vclume  element  5V.  It  is 


5N  =  S(z )  T(z  ,l)  -5V  (21) 

'4- 

where  ,  is  the  solid  angle  for  emission  which  enters  the  detector.  S  is 
assumed  to  have  an  isotropic  behavior  which  leads  to  the  presence  of  1/4-. 

We  wish  to  integrate  ?N  over  A^  since  this  is  the  actual  observing  cross 
sectional  area  at  distance  R.  We  assume  S  to  be  constant  over  A_.  Letting 
f N ’  refer  to  the  integrated  quantity,  it  is 

fN '  =  S(z)  T ( z  , „ )  A  * s  iv_  (22) 

4- 


We  now  wish  to  modify  fN'  to  remove  its  dependence  on  any  particular  detector. 
To  do  so,  we  divide  6N '  by  A^f.^  where  Ad  is  the  area  of  the  detector  and 
is  the  detector's  solid  angle  field  of  view.  Since  the  modified  fN'  has  the 
units  of  an  incremental  intensity,  let  it  be  denoted  by  51.  It  is 


51 


A 

E 


T (z  ,u)6s 


(23) 
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The  quantity  ./  ^  A  /A^  is  unity  as  can  be  seen  by  noting  that  .  =  A^/P 
and  d  =  A  /R^.  Replacing  ds  by  dz/u.  Equation  (23)  becomes 


■:I  =  ^  T  (z,u  )  dz/u 


(24) 


which  can  be  eauivantly  expressed  as 


'I  =  ^  T (  t / m )  dt/p 


(25) 


Integration  of  51  gives  I  as 


shown  in  Equation  (15). 


6.4  Forms  of  SQ 

Expressions  will  follow  for  the  initial  volume  emission  rate 
S  (z).  In  the  ionosphere,  we  are  here  interested  in  excited  state  pro¬ 
duction  by  electron  impact.  Let  us  introduce  P,  the  volume  excitation 
rate  which  equals  SQ(z)  only  if  quenching  is  unimportant.  The  expression 
for  P  due  to  electron  impact  excitation  is 


P£k(z)  =  r  n  ,(z)  J  o,.k(E)  4>  (z,E)  dE 


(26) 


3 

with  units  of  excitations/cm  -s.  The  sum  is  over  species  and  applies  to 

4 

certain  atomic  states  such  as  N(  P)  which  can  be  produced  by  either 
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or 


N(4S)  +  e  -«•  N(4P)  +  e 
N2(X)  +  e  -  N(4P)  +  N  +  e. 

States  of  interest  to  us  will  be  effectively  produced  by  electron  impact 
on  a  single  type  of  neutral  particle.  The  terms  k  and  4-  in  Equation  (26) 
are  the  electron  impact  cross  section  and  spherical  flux  for  either  auroral 
or  photoelectrons.  The  index  k  refers  to  the  particular  state  which  may 
be  an  atomic  state,  an  electronic  state  of  a  molecule,  or  a  vibrational 
state  of  a  molecule. 

Let  us  now  relate  SQ  to  P^  and  in  the  process  we  shall  delete 
the  indices.  For  an  atomic  state  unaffected  by  quenching, 

SQ(z)=P(z).  (27) 


With  quenching,  the  relationship  is 


where  A  is  the  transition  rate  (s'1)  for  the  emission  of  interest,  AT  is  the 

3-1 

total  rate  from  the  upper  state,  and  q.  is  the  quenching  coefficient  (cm  -s  ) 
+  h 

for  the  v  species. 


For  molecules,  we  are  interested  in  emission  within  a  band.  Let  us 

assume  P  refers  to  the  excitation  rate  for  a  particular  v'  level  of  the  excited 

electronic  state.  Excitation  may  be  a  combination  of  direct  excitation  from 

3 

the  ground  state  and  cascading.  In  the  case  of  the  A  i  state  of  N£,  cascading 
is  the  dominant  excitation  mechanism  (see  Cartwright,  1978).  In  the  absence 
of  quenching,  S0  is  related  to  Py.  by 
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so(z)  'rnrh  V(2) 

ii  v  v 


(2-6) 


where  Av,v„  is  the  transition  rate  from  level  v’  of  the  uprer  state  to  level 
v"  of  the  lower  state.  With  quenching,  the  relationship  is 


So(z)  = 


Ay'y"  Pyl 


Av'v"  +  1  Vi<z) 


(30) 


similar  to  Equation  (28)  for  atomic  emission. 


For  Rayleigh  scattering,  we  are  interested  in  two  different  ex¬ 
pressions  for  SQ  -  one  for  the  unidirectional  solar  source  and  the  other 
for  diffuse  auroral  emissions.  For  the  solar  case,  SQ  has  the  form 


s0(z)  -  *F0(>)e 


_  I  ^  X 

v  lR  ' 


+  / 
t- 1 1 , 


ncr 


(31) 


where  rF  is  the  unattenuated  solar  flux  at  some  given  wavelength  '  in  photons/ 

2  0 

cm  -s-A,  n  is  the  total  neutral  scattering  density,  and  is  the  Rayleigh  cross 
section.  By  normalizing  SQ  to  a  unit  solar  flux,  Sq(t)  is  simply 


So(tR}  = 


•(td  +  t)/uc 


(32) 


which  is  a  convenient  form  for  comparisons  of  S0  and  S  from  one  wavelength 
to  another. 
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The  form  of  SQ(z)  for  an  auroral  radiation  source  is 


SQU)  =  noR 


Kfi)  e 


(Tr  +  t)/p 

K  da 


'3 


where  I  is  the  auroral  intensity  for  some  feature  of  interest  as  it  leaves 
the  ionosphere  in  the  downward  direction  given  by  fi.  Here,  z,  as  in  Equation 
(31).  refers  to  altitudes  in  the  lower  atmosphere.  We  may  think  of  the 
integration  region  in  Equation  (33)  as  being  over  2tt  with  1(a)  properly  set 
to  zero  for  directions  from  which  no  auroral  luminosity  is  emitted.  In 
reality,  there  can  be  a  variety  of  directional  luminosity  distributions 
from  those  giving  limited  ranges  such  as  for  a  single  arc  to  more  hemispherical 
distributions  as  can  be  found  under  weak,  diffuse  conditions. 


6.5  Parameter  Values 

The  basic  parameters,  already  Introduced,  are 

2 

1)  oR(A),  the  Rayleigh  cross  section  in  cm 

2)  ir(A),  the  corresponding  optical  depth 

3)  Op(A),  the  pure  absorption  cross  section 

4)  tp( A ) ,  the  corresponding  optical  depth 
and 

2 

5)  ttF  (A),  the  unattenuated  solar  flux  in  photons/cm  -s-A. 
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,  we  use  the  expression 


For  lR 


Op  ( A ) 


3.9  x  10~28 

.  a 
A 


cm 


3 


with  a  given  by 


a  =  3.916  +  .074  A  +  ~ 

(Frolich  and  Scaw,  1980).  The  units  of  A  are  10‘4  A.  Figure  7  displays  -R 
between  1800  and  3200  A. 

Pure  absorption  is  caused  by  Op  and  0^  over  the  wavelength  range 
of  interest.  For  0^,  -p  is  shown  in  Figure  8  from  1700  to  ^  2600  A,  taken 
from  data  gathered  by  Sullivan  and  Holland  (1966).  The  region  shown  below 
2000  A  is  dominated  by  Schumann-Runge  band  absorption.  Above  2000  A, 
absorption  is  associated  with  the  Herzberg  continuum.  Figure  9  shows  ;p 
for  Op  between  2000  and  3000  A,  again  taken  from  Sullivan  and  Holland 
(1966).  The  dominant  absorption  in  this  region  is  Hartley  continuum 
absorption . 


The  depths  tr  and  tp  at  z  =  0  km  are  shown  in  Figure  10.  Altitude 
profiles  at  selected  wavelengths  appear  in  Figure  11.  We  note  from  Figure  10 

that  t  becomes  small  as  well  as  less  than  longward  of  n-  3000  A.  As  will 
oe  seen,  this  leads  to  a  noticeable  increase  in  Rayleigh  scattering  long- 
ward  of  this  wavelength.' 
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'SS  SECTION 
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An  example  of  a  solar  spectrum  is  shown  in  Figure  12  from  1800  A 
to  3200  A,  (Banks  and  Kockarts,  1973A).  The  only  point  we  wish  to  make  here 
is  that  even  weak  Rayleigh  scattering  at  wavelengths  such  as  3000  A  can 
swamp  i onospherical ly  produced  intensities  as  observed  from  satellite 
altitudes  in  the  daytime  due  to  the  strength  of  the  solar  flux.  More 
quantitative  information  will  follow  in  the  next  subsection. 


6.6  Results 

We  begin  by  showing  the  effect  of  multiple  scattering  over  the 
wavelength  range  2000  -  4500  A.  This  is  done  in  Figures  13  and  14  which  show 
S  ( t )  and  S(')  at  altitudes  of  96  and  47  km  (D.E.  Anderson,  private  communica¬ 
tion).  S( : )  was  obtained  from  Equation  (20)  with  S0(-r)  given  by  Equation  (32). 

SQ  thus  refers  to  the  initial  scattering  of  solar  radiation  (with  ^  set  to 
unity).  A  similar  -  dependence  is  exhibited,  however,  for  the  initial  Rayleigh 
scattering  of  auroral  radiation  and  thus  this  discussion  applies  to  either 
situation.  By  comparing  S  to  SQ  within  either  Figure  13  or  14,  we  see  that  multiple 
scattering  only  becomes  important  above  3000  A.  Below  this  wavelenqt1’,  pure 
absorption  is  strong  enough  by  0£  and  0^  to  prevent  radiation  *rom  effectively 
penetrating  to  altitudes  where  the  condition  r  :  1  is  met. 

We  now  present  results  for  the  backscatter  of  solar  radiation. 

We  use  Equation  (15)  to  calculate  4-1,  the  apparent  column  emission  rate. 

The  T  function  is  just  e-^  Then,  setting  S  to  SQ  as  aiven  by  Equation 

(32),  4rl  is 

4- 1 R  (;. ,  „s )  =  rFofs  ph/cm2-s-A  (36 


45 


S(i)  (ph./cm  -unit 


with  f  given  by 


o 


Note  that  the  Rayleigh  phase  function  dees  not  appear  because  of  'he  iso¬ 
tropic  scattering  approximation  justified  earlier  in  this  discussion. 

Figure  15  shows  f  and  4~IR  as  functions  of  wavelength  from  14'JO 
to  3200  A.  For  these  results,  p  =  =1.  The  peak  in  f&  at  200C  A  is  due 

to  a  minimum  in  the  combined  pure  absorption  by  Op  and  0^.  It  should  be 
noted  that  structure  in  f  will  be  present  below  2000  A  due  to  structure  in 
Schumann-Runge  band  absorption  which  is  responsible  for  the  exhibited  de¬ 
crease.  The  decrease  in  f^  longward  of  2000  A  is  caused  by  Op  absorption. 
The  dramatic  rise  near  3000  A  results  from  this  absorption  decreasing  to 
the  point  that  radiation  can  penetrate  deeply  into  the  atmosphere.  We 
have  included  multiple  scattering  effects  above  3000  A  in  an  approximate 
way.  This  has  been  done  by  scaling  4ml  as  calculated  in  Equation  (36) 
by  S ( t )/S  h)  with  t  given  near  its  mid-point  value. 

The  profile  for  4~IR  was  obtained  using  ttFo  from  Figure  12 
We  observe  that  more  than  a  ki 1 o- Rayl ei gh  per  Angstrom  is  backsca ttered 
above  1900  A.  Thus,  satellite  observations  must  be  restricted  to  wave¬ 
lengths  shorter  than  this  since  ionospherically  produced  emission  features 
above  1900  A  will  be  swamped  by  the  given  Rayleigh  intensity  component. 

Our  last  results  are  for  the  Rayleigh  component  of  auroral 
observed  intensities.  As  we  shall  see,  the  observational  restrictions 
with  regard  to  wavelength  are  less  severe  here  since  there  is  no  strong 
external  radiation  source  as  there  is  in  the  daytime.  Again  we  replace 
S  by  SQ  with  SQ  here  given  by  Equation  (33).  To  obtain  an  upper  limit 
to  SQ,  we  assume  the  auroral  luminosity  I  (ft)  to  be  uniform  and  independent 
of  P  across  the  sky.  Then  Equation  (33)  becomes 
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Figure  15  Functions  f  and  4 
of  Solar  Radiation 


5 


Sq(t)  =  2nl  j  E2  (t  +  t) 


(38) 


where  E^  is  the  second  exponential  integral.  Inserting  this  into  Equation 
(15) ,  4tIr  becomes 


4tIr  (  j )  =  4-rrIj  fft 


ph/cm-s 


(39) 


with  f^  given  by 


J  °  E 2 ( T  +  t)  e'(T  +  t)/u  dt/u 


(40) 


The  total  column  emission  rate  4ttI  is  then 


4ttI(p)  =  4 rrl j  (1  +  fA) 


(41) 


Figure  16  shows  from  1800  A  to  3200  A.  The  behavior  exhibited  is  similar 
to  that  of  f$  and  as  with  that  function,  multiple  scattering  effects  have 
been  accounted  for  at  the  longer  wavelengths.  We  see  that  Rayleigh  scattering 
is  unimportant  below  n,  3200  A. 
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Section  7 


Predictions  for  Comparison  with  Data  from  the 
1981  AFGL  Auroral  £  Rocket  Program  and 
from  the  ISIS  II  Satellite 


In  this  section,  we  address  two  unrelated  experimental  programs 
which  have  provided  various  types  of  data  under  auroral  conditions.  The  first 
to  be  considered  will  be  the  AFGL  Auroral  E  Rocket  Program  which  was  carried 
out  on  March  7,  1981.  Four  rockets  were  launched  from  the  Poker  Fiat  Research 
Range  near  Fairbanks,  Alaska  into  the  continuous  auroral.  Observations  were 
made  of  particles,  fields,  and  optical  emissions.  The  second  program  involves 
the  ISIS  II  satellite  and  of  specific  interest  here  are  observations  which  were 
made  with  soft  particle  spectrometers  (see  Klumpar  and  Heikkila,  1982). 


7.1  Predictions  for  the  Continuous  Aurora 

The  quantities  of  interest  here  are  the  following: 

1)  primary  and  secondary  electron  spectra 

2)  ion  densities 

and  3)  selected  optical  intensities  (01  1356  A,  01  5577  A,  01  2972  A, 

N*  IN  3914  A,  N2  2P  3371  A,  several  LBH  bands,  and  finally 
the  N2  Vegard-Kaplan  (VK)  band  at  2762  A). 

They  have  been  calculated  using  a  Jacchia  (1977)  model  atmosphere 

with  the  incident  electron  spectrum  represented  by  a  1  keV  Maxwellian  distribu- 

o 

tion.  The  energy  content  Q  is  1.5  ergs/cm  -s  and  isotropy  is  assumed  over  the 
downward  hemisphere.  The  chosen  representation  is  close  so  that  deduced  from 
coincident  Chatanika  radar  data  (Robinson  and  Vondrak,  1981)  and  from  AFGL  Air¬ 
borne  Laboratory  data  (j.A.  Whalen,  private  communication,  1981).  Primary  electron 
data  were  not  available  during  the  study  for  direct  input  into  the  transport 
calculations. 
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We  begin  with  the  electron  transport  results  in  Figure  17.  Shown 
are  spherical  fluxes  (e/cm^-s-eV-4irsr)  at  the  altitudes  300  km  (upper  boundary), 

135  km,  and  120  km  over  the  energy  range  1.5  eV  to  5  keV.  w»  note  that  we  h*»e 
allowed  the  incident  spectrum  to  depart  from  a  Maxwellian  below  ^  200  eV.  if 
starts  with  a  true  Maxwellian  it  will  be  seen  to  evolve  to  the  spectrum  shown  at 
300  km  due  to  secondary  electron  production.  Th-;s  occurs  before  any  noticeable 
energy  degradation  of  the  primary  spectrum  has  taken  place  above  a  few  hundred  eV. 

VJe  thus  choose  to  start  with  this  slightly  evolved  version  of  the  incoming  spectrum. 

The  minima  occuring  at  %2.5  eV  are  due  to  vibrational  excitation. 

In  this  low  energy  region,  some  of  the  loss  is  also  caused  by  plasma  electrons  as 
well  as  transitions  to  excited  states  with  low  energy  thresholds.  The  plasma 
density  needed  for  describing  loss  to  the  plasma  was  chosen  to  be  consistent  with 
the  incident  electron  spectrum. 

The  results  shown  in  Figure  17  are  based  on  a  particle-particle  inter¬ 
action  description.  They  show  stronger  variations  below  100  eV  than  have  some¬ 
times  been  observed  (see,  e.g.  Sharp  and  Hays,  1972  and  Feldman  and  Doering,  1975). 
Such  discrepancies  have  been  suggested  to  arise  from  wave-particle  interactions 
not  accounted  for  in  the  calculations  (see  Papadopoulis  and  Coffey,  1974).  Our 
point  here  is  to  acknowledge  a  potentially  important  role  of  plasma  turbulence  and 
note  under  what  conditions  our  own  calculations  have  been  made.  It  seems  un¬ 
likely,  however,  that  the  calculated  secondary  spectrum  is  in  serious  error  when 
taken  over  the  entire  altitude  range  of  interest.  This  is  based  on  its  associat¬ 
ed  emission  rates  and  their  agreement  with  observations  (see  Meier  et  al.,  19821. 

Figure  18  shows  calculated  ion  densities  for  N0+,  02+,  as  well  as 

their  sum  which  is  equated  to  the  electron  density  n  .  These  results  come  from 
our  time  dependent  chemistry  model  which  does  not  include  the  effect  of  ion 
transport.  The  calculation  was  continued  until  chemical  equilibrium  was  achieved. 
The  behavior  of  the  N0+  and  0^  densities  are  sensitive  to  the  amount  of  NO  pre¬ 
sent.  The  chosen  NO  density  altitude  profile  may  be  found  in  Table  1  of  the 
appendix. 
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Energy  (eV) 

Calculated  auroral  electron  spectra.  The  incident 
spectrum  is  that  shown  at  300  km  which  has  the  dependence 
of  a  1  keV  Maxwellian  distribution  above  n.  200  eV  where 
most  of  the  energy  resides. 


The  densities  in  Figure  18  nay  be  compared  with  various  data.  We 
observe  the  same  basic  behavior  in  ng  deduced  by  Robinson  and  Vondrak  (1981) 
although  their  profile  is  'u  50%  higher.  This  is  not  a  serious  difference  given 
that  we  do  not  have' a  direct  measure  of  what  the  actual  energy  content  of  the 
incident  electron  spectrum  was.  Figure  IS  provides  a  comparison  between  our  re¬ 
sults  and  the  data  of  Swider  and  Narcisi  (1981)  for  fractional  abundances  of  ions. 
The  data  were  obtained  with  an  ion  mass  spectrometer.  The  overall  agreement  is 
good  with  the  exception  of  for  which  our  results  are  ^  3  times  larger.  We 
have  not  yet  investigated  the  cause  of  this  discrepancy. 

We  now  consider  the  optical  emission  features  listed  earlier.  Column 
emission  rates  (to  be  referred  to  as  intensities)  in  kilo- Ray 1 e i ghs  (kR)  will  be 
presented  for  viewing  in  the  zenith  direction.  Data  for  the  chosen  features  were 
obtained  for  this  viewing  direction  (R.E.  Huffman,  private  communication,  1981, 
R.A.  Van  Tassel,  private  communication,  1981).  We  begin  with  the  features  IN 
3914  A  and  01  5577  A  whose  intensities  are  among  the  most  prominent  in  auroras. 
Figure  4  shows  altitude  profiles  of  these  intensities  (previously  presented  in 
our  discussion  of  0(1S)  chemistry).  Both  features  yield  intensities  in  excess  of 
1  kR  for  the  given  aurora  with  being  o,  50%  larger  at  low  altitudes.  This 

is  consistent  with  numerous  observations  which  typically  show  the  ratio  4^1.^/ 

4’t 1 3914  to  be  between  1  and  2. 

Figure  20  shows  intensities  for  01  2972  A,  2P  3371  A,  and  VK 
2762  A.  The  differences  between  the  2762  A  and  3371  A  profiles  are  caused  by 
quenching  which  affects  the  A  state  leading  to  2762  A  emission  but  not  the 
C  state  which  leads  to  3371  A  emission.  Intensities  for  selected  LBH 
bands  are  shown  in  Figure  21.  The  profile  labeled  1493  A  contains  contributions 
from  NI  1493  A  and  the  (  3,3)  band  at  the  same  wavelength.  The  other  profiles 
contain  contributions  from  pairs  of  bands  as  indicated  in  the  figure.  Pure 
absorption  by  0^  is  responsible  for  the  decreases  at  low  altitudes.  It  is 
weaker  near  1425  A  than  at,  say,  1384  A  which  leads  to  the  smaller  decrease 
below  120  km  for  the  combined  1426  A  -  1430  A  intensity  profile.  The  LBH 
1354  A  intensity  is  shown  in  Figure  22  along  with  that  for  01  1356  A  as  well  as 
their  sum.  For  the  given  incident  electron  spectrum,  we  observe  that  the  LBH 
contribution  to  the  summed  intensity  is  minor.  It  will  become  more  important 
as  the  hardness  of  the  electron  spectrum  increases  since  in  this  situation,  more 
of  the  energy  becomes  deposited  in  N^  relative  to  0. 
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Figure  20  Calculated  zenith  viewing  intensities  for 
conditions  described  in  text. 


Fiaure  21  Calculated  zenith  viewing  intensities  for 
conditions  described  in  text. 
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Figure  22  Calculated  zenith  viewing  intensities  for 
conditions  described  in  text. 
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7.2  Election  Backscatter  Prediction  for  ISIS  II  Electron  Spectrometer 
Data  Analysis 

Two  sets  of  auroral  electron  spectral  data  were  receivv  :>c  >! 

with  the  request  to  calculate  backscattered  spectra.  The  interest  in  these 
particular  data  come  from  their  beam  like  quality  in  the  downward  herr.ispnere 
(strongly  field  aligned  at  certain  energies)  and  large  values  at  low  energies 
over  the  upward  hemisphere.  The  behavior  of  some  of  the  data  strongly  suggests 
the  existence  of  parallel  electric  fields  below  the  satellite.  Klumpar  and 
Heikkila  (1982)  have  specifically  addressed  this  issue  using  data  from  the  ex¬ 
periment  we  are  considering  here.  Our  calculations  were  requested  to  determine 
the  extent  of  backscatter  due  to  collisions  so  as  to  better  quantify  the  inter¬ 
action  of  electrons  (plasma  and  energetic)  with  electric  fields  beneath  the 
satel 1 i te . 

Figure  23  shows  an  example  of  the  data  used  to  specify  the  incident 
electron  spectrum.  The  left  and  right  panels  contain  spectra  at  pitch  angles 
of  V  and  14'  respectively  and  illustrate  the  beam  like  nature  of  the  data  be¬ 
low  a  few  keV.  The  altitude  at  which  the  observations  were  made  is  1400  km. 
Transport  calculations  were  done  for  two  sets  of  data  which  we  shall  designate 
by  the  labe1s  A  and  5.  For  case  A,  data  are  available  at  pitch  angles 
of  7  .  14',  27',  48',  68',  and  88°.  Case  B  data  are  available  at  11°.  31c  , 
and  33  .  (it  should  be  noted  that  data  do  exist  at  angles  above  90  but  were 
rot  received  by  us  for  this  analysis.)  The  pitch  angle  is  here  referenced  to 
the  downward  direction  along  the  geomagnetic  field  lines.  The  data  a  53'  for 
case  B  are  questionable  since  numerous  data  dropouts  occur  and  are  otherwise  of 
low  magnitude  with  considerable  fluctuation.  Case  A  results,  therefore,  deserve 
the  greater  emphasis. 

Figure  24  shows  our  smoothed  representation  of  the  data  and  the 
calculated  backscattered  spectra  at  selected  energies.  The  calculations  were 
performed  with  many  more  energy  points  than  shown  in  the  figure.  Twenty  non- 
umformly  spaced  \i  values  (u  is  the  cosine  of  the  pitch  angle)  were  used  for 
both  cases.  At  high  energies,  the  backscattered  spectrum  is  weak  as  expected  sine 
fo-ward  scattering  is  strong  and  few  electrons  are  available  at  higher  energies 
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Figure  24  Incident  (u  >  0)  and  calculated  backscattered  (u  <  0) 

electron  spectra  based  on  ISIS  II  data  such  as  appear  in  the 
previous  figure. 


which  might  contribute  after  considerable  energy  degradation.  At  low  energies, 
the  upward  flux  is  appreciable  due  mainly  to  the  rate  at  which  secondary  electrons 
are  produced.  If  one  examines  the  data  in  the  paper  of  Klumpar  and  Heikkila  (1982), 
however,  it  is  clear  that  the  amount  of  backscatter  exhibited  in  our  results  is 
insufficient  at  times  to  explain  observed  low  energy  upward  electron  fluxes,  thus 
suggesting  some  acceleration  mechanism  in  the  upper  ionosphere. 


.St ‘Oti^K  8 

UV  Emissions  and  the  Electron  Density  in  the  Daytime  Ionosphere 


8. 1  Background 


The  problem  we  address  in  this  section  has  to  do  with  properties 
of  UV  emissions  and  the  electron  density  in  the  daytime  ionosphere.  These  quant¬ 
ities  vary  with  geographical  location,  local  time,  and  from  day  to  day.  The 
causes  of  variation  in  emissions  other  than  local  time  variation  are  1)  changes 
in  the  solar  flux  producing  ionization  and  2)  changes  in  the  neutral  atmosphere 
above  '  100  km.  These  changes  also  directly  affect  the  electron  density.  The 
presence  of  plasma  turbulence,  winds,  and  large  scale  weak  electric  fields  will 
also  affect  the  electron  density,  particul ari ly  in  the  F  region.  Here,  we  will 
restrict  ourselves  to  emission  and  electron  density  variations  caused  by  composi¬ 
tion  and  solar  flux  changes. 

We  are  particularly  interested  in  knowing  how  useful  UV  emissions 
are  in  determining  the  daytime  electron  density  profile  in  the  E  and  F  regions. 

We  leave  to  a  later  time  considerations  of  how  accurately  such  a  determination  must 
be  for  the  technique  to  be  judged  useful.  To  begin  addressing  the  problem,  code 
PEGFAC  has  been  placed  on  the  AFGL  computing  system  and  coupled  to  code  CHEM.  We 
are  now  beginning  to  obtain  results  using  these  codes.  These  results  include 
photoel ectron  fluxes,  photoelectron  produced  excitation  and  ionization  rates, 
primary  ionization  rates  caused  directly  by  the  solar  flux,  densities  of  all 
major  chemical  species,  and  optical  intensities. 

PEGFAC  refers  to  jjhotoel ectron  g- factors  which  are  the  parameters 
that  provided  the  initial  motivation  for  us  to  develop  the  code.  A  g-factor  is 
defined  to  be  some  volume  excitation  or  ionization  rate  divided  by  the  density  of 
the  parent  species  from  which  the  excited  species  originated.  For  photoelectron 
produced  excitation,  it  is  als"  the  integral  over  energy  of  the  excitation  cross 
section  times  the  photoelec tron  flux.  Code  PEGFAC  provides  for  the  calculation 
cf  this  flux  for  various  solar  conditions  and  model  atmospheres  as  well  as 
calculation  of  a  larqe  number  of  excitation  rates  and  their  corresDondi no 
g-factors. 
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Code  CHEM  is  our  auroral  chemistry  code  which  calculates  the 
densities  of  the  important  chemical  species  as  well  as  numerous  intensities 
of  molecular  band  systems  and  atomic  lines.  The  code  now  allows  for  diffusion 
of  0+  and  thus  enables  us  to  model  the  electron  density  into  the  F  region. 

Some  code  modifications  have  been  necessary  to  adapt  CHEM  to  the  daytime  problem 
which  are  now  complete. 

In  the  following  subsections,  some  information  will  be  given 
concerning  the  model  incorporated  into  PEGFAC.  Calculated  intensities  and 
electron  densities  will  also  be  given  as  we  begin  to  address  the  problem  as 
stated  at  the  beginning  of  this  subsection. 

As  an  historical  note,  a  significant  part  of  the  coding  leading 
to  PEGFAC  was  done  under  DARPA/NRL  sponsorship  to  address  the  problem  of  secondary 
electron  energy  deposition  during  and  following  the  passage  of  high  current-high 
energy  electron  beams  in  air.  Unsponsored  work  by  one  of  us  (D.J.S.) 
then  led  to  the  initial  untested  version  of  PEGFAC.  Most  of  this  effort  was 
directed  to  development  of  a  photoelectron  source  code.  Valuable  assistance  was 
given  by  R.R.  Meier  and  colleagues  at  NRL  in  providing  the  necessary  detailed 
information  on  photoionization  cross  sections  and  solar  fluxes.  Testing  and 
initial  application  of  PEGFAC  was  then  carried  out  under  NASA  sponsorship  through 
the  Aerospace  Corporation.  Part  of  the  descriptive  information  to  follow  is 
taken  for  the  documentation  of  this  effort.  The  results,  however,  come  from  this 
work . 


8 . 2  Computational  Model 

The  photoelectron  flux  is  calculated  from  an  integral  equation 
in  the  local  approximation.  The  equation  is  basically  that  crnsidered  by 
Strickland  et  al.  (1976)  for  auroral  electrons  and  by  Oran  and  Strickland  (1978) 
for  photoelectrons  without  the  ud-f/dz  term.  Deleting  the  transport  term  greatly 
simplifies  the  calculations  and  is  valid  for  our  purposes  over  the  altitude  range 
of  important  impact  excitation.  We  have  compared  our  solution  with  one  generated 
using  the  Onan-Strickland  code  and  observe  dose  agreement  to  300  km. 
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The  integral  equation  is 


__  __  __  r  max 

0  *  -  n4(2)ap_(E  )0(z,E)  +  /  ‘\k(E‘  ,E);(z,E'  )d£ 

fi.  ^  k  ^ 

+  np  Ie  CLp(E)<t(z,E)]  +  S{z,E)  (42) 

with  terms  defined  as  follows 

2 

4  photoelectron  flux  in  e/cm  -s-eV 

3 

5  source  spectrum  in  e/cm  -s-eV 

n,  density  of  £th  neutral  species  (N^.O^  and  0  treated) 

f-  h  O 

c„  total  inelastic  cross  section  of  species  in  cnr 

t  h  th 

o. ,  inelastic  cross  section  for  kLn  process  involving 

■  k  2 

species  in  cm  /eV 

np  plasma  density 

Lp  loss  function  for  energy  loss  to  plasma- 


We  are  interested  in  a  variety  of  volume  excitation  rates  (the 
term  "excitation"  will  be  used  at  times  for  both  excitation  and  ionization)  for 
providing  chemistry  source  terms  and  volume  emission  rates.  These  excitation 
rates  come  from  photoelectron  impact  excitation  and  in  some  cases  from  direct 
solar  photon  excitations.  The  rates  of  concern  here  produced  by  both  processes 
are  ionization  rates.  The  direct  solar  part  is  given  by 
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(43) 


P$(z,u0)  =  n(z)  g$(z,uo)  cm^-s'1 


with  u  being  the  cosine  of  the  solar  zenith  angle  and  the  solar  ionizatior 
g-factor  given  by 


/-t, (z)/un 

7Fxe  d. 


The  terms  mF^,  t^,  and  refer  to  the  unattenuated  solar  flux  (cm  -s-A)  ,  the 
total  photo-absorption  optical  depth,  and  the  photo-ionization  cross  section 
for  some  specific  ion  state. 

Since  P$  refers  to  some  specific  ion  state,  it  is  appropriate  to 
include  various  subscripts  in  Equations  (43)  and  (44).  We  have  chosen  to  leave 
them  off  for  the  sake  of  simplicity  in  notation.  It  should  be  noted,  however, 
that  a  sum  will  be  involved  in  Equation  (43)  if  more  than  one  species  leads  to 
production  of  the  given  ion  species. 


The  volume  excitation  rate  by  photoelectrons  for  some  specific 


process  is 


PpE(z,u0)  =  n(z)  gpE(z,u0)  cm"3-s_1 


with  the  photoelectron  g-factor  given  by 


/max 

(z,E)r>(E)dE 

o 
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where  W  is  the  excitation  threshold  energy  and  subscripts  uQ  has  been  added  to 
c  to  designate  that  4  has  this  dependence.  Again,  for  simplicity  in  notation, 
we  chose  not  to  show  this  dependence  in  Equations  (42)  and  (43). 

The  next  quantity  to  be  defined  is  the  intensity.  We  shall  con¬ 
sider  its  computational  form  for  a  plane-paral lei  geometry.  Designated  as 
I(z.u),  it  is 


I(z,u)  =  (4-r) 


7 


S(z 1 )  T(z,z' ,u)dz '/ 


ph/cm  -s -sr 


(47) 


where  S(z)  is  the  total  volume  emission  rate  and  T  is  the  transmission  function 
between  altitudes  z  and  z'  along  a  slant  path  specified  by  p,  the  cosine  of  the 
look  direction.  We  have  not  shown  the  integration  limits  since  they  depend  on 
whether  j  is  positive  (the  convention  here  is  that  u  >  0  corresponds  to  nadir 
viewing)  or  negative  (looking  up).  The  emission  rate  includes  a  multiple  scattered 
contribution  if  the  medium  is  optically  thick.  Here,  we  will  consider  only  media 
which  are  sufficiently  thin  that  S(z)  can  be  adequately  represented  by  S  (z),  the 
emission  rate  with  no  multiple  scattering  contribution. 

The  final  quantities  to  be  defined  are  the  densities  determined 
from  chemistry  modeling.  These  include  densities  for  both  neutral  and  ion 
species  in  ground  and  excited  states.  The  applied  rate  equation  for  each  species 
is 


(48  1 


where  p ■  and  are 
shall  be  interested 


t  h 

the  production  and  loss  rates  for  the  v  species.  We 
in  the  steady  state  solutions  to  Equation  (48).  To  solve 
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the  coupled  rate  equations,  we  have  applied  code  CHEM,  our  auroral  chemistry 
code.  Transport  of  0  is  treated  in  CHEM  and  consequently  we  may  examine 
ionospheric  properties  in  both  the  E  and  F  regions. 


8.3  Description  of  Input/Output  Information. 

The  calculation  of  photoelectron  spectra  and  associated  volume 
excitation  rates  requires  a  large  body  of  input  data.  In  this  subsection, 
we  will  note  the  extent  of  the  needed  input  parameters. 


The  input  information  may  be  catagorized  as  follows 

•  model  atmosphere 

•  solar  spectrum  (EUV  to  X-ray) 

•  photoabsorption  cross  sections 

•  partial  photoionization  cross  sections 

•  inelastic  cross  sections  for  modeling  the 
photoelectron  energy  degradation 

•  loss  function  for  energy  loss  to  the  plasma 

and  •  excitation  cross  sections  for  specifying 
excitation  processes  of  interest. 


Some  of  these  quantities  have  already  been  discussed.  Never¬ 
theless,  we  will  briefly  address  each  of  them  here  noting  either  their  source 
or  extent.  The  model  atmosphere  comes  from  Jacchia  (1971)  and  is  generated  in 
the  code  once  the  exospheric  temperature,  Tw,  is  specified.  We  allow  for  a 
scaling  of  the  0  density  through  a  factor  appearing  in  the  input  data. 

The  solar  flux  values  come  from  Torr  et  al.  (1979)  and  Donnelly 
and  Pope  (1973).  In  a  given  run,  one  of  five  available  spectra  is  selected 
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depending  on  the  degree  of  solar  activity  to  be  modeled.  A  combination  of 
line  and  continuum  band  fluxes  totaling  38  in  number  is  provided  with  the 
Torr  et  al.  spectra.  We  have  replaced  their  short  wavelenqth  band  fluxes 
with  line  fluxes  from  Donnelly  and  Pope  (1973).  A  scaling  factor  has  been 
introduced  to  allow  experimentation  on  our  part  for  lack  of  knowledge  on  how 
the  soft  X-ray  fluxes  vary  with  time. 

The  total  absorption  cross  sections  are  used  to  specify  the 
attenuated  solar  spectrum.  They  come  from  Torr  et  al.  (1979).  The  partial 
ionization  cross  sections  come  from  Kirby  et  al.  (1979)  and  selected  papers 
referenced  therein.  Photoionization  is  modeled  for  five  states  of  N?+,  ten 

■f  +  ^ 

states  of  C>2  ,  and  five  states  of  0  .  The  ionization  thresholds  span  an  energy 
range  from  12.1  to  25  eV. 

There  are  two  sets  of  electron  impact  cross  sections.  The 
first  set  defines  the  energy  loss  matrix  elements  while  the  second  set  con¬ 
tains  members  for  excitation  to  states  of  specific  interest  to  dayglow  studies. 

In  terms  of  energy  loss,  the  latter  is  a  subset  of  the  first.  Most  of  the 
applied  cross  sections  come  from  the  previous  work  of  Strickland.  These  may 
be  found  in  the  papers  by  Strickland  et  al .  (1976)  and  Oran  and  Strickland 

(  1 9 7 8 ) .  Recent  members  have  been  ,  added  for  processes  such  as  dissociative 
ionization  of  09  leading  to  011  834  A.  Most  of  these  have  come  from  measurements 
by  E.C.  Zipf  and  colleagues. 

The  final  parameter  on  the  above  list  is  the  loss  function  for 
energy  loss  to  the  plasma.  It  should  be  noted  that  such  loss  is  not  important 
to  the  study  of  the  UV  dayglow.  Electron  impact  excitation  leading  to  UV 
emission  is  primarily  at  electron  energies  above  10  eV  while  plasma  absorption 
of  the  photoelectron  energy  which  affects  the  flux  spectrum  occurs  below  this 
energy.  In  spite  of  this,  we  include  the  effect  of  plasma  energy  loss  so  as 
not  to  exclude  the  low  energy  region  from  future  studies.  The  applied  loss 
function  comes  from  Schunk  and  Hays  (1971). 
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We  will  now  briefly  describe  the  output  information.  The  flux 
:(z,E)  is  currently  calculated  for  n.  40  altitudes  between  400  and  100  km  and 
for  as  many  as  75  energies.  Approximately  40  excitation  rates  and  45  g-factors 
are  specified  on  the  altitude  grid.  These  refer  to  ionization,  vibrational 
excitation,  electronic  state  excitation  of  the  states  producing  important  energy 
loss,  and  additional  processes  leading  to  UV  emission  at  selected  wavelengths. 

The  excited  state  species  producing  the  emission's  include  N,,,  N^+,  0,  0  ,  and  N. 


8 . 4  Calculated  Photoelectron  Flux  and  Volume  Excitation  Rates 

An  example  of  the  calculated  photoelectron  flux  as  obtained  from 
Equation  (42)  is  shown  in  Fiqure  (25)  together  with  data  by  Lee  et  al.  (1980).  Both 
spectra  apply  to  low  solar  activity  and  are  seen  to  be  in  good  agreement  with  one 
another.  The  applied  solar  spectrum  is  shown  in  Figure  26  with  conditions  given 
in  the  caption. 

Figure  27  shows  exam^cs  of  photoelectron  produced  volume  excitation 
rates  for  (a  ^g)  and  0  (5S).  Conditions  are  the  same  as  applied  to  the  in¬ 

formation  in  Figures  25  and  26.  Solar  zenith  anqles  are  considered  over  the  ranae 
from  0'  to  90'.  The  given  states  were  chosen  since  we  are  currently  examining 
their  associated  emissions  which  occur  in  the  LBH  bands  and  at  1356  A. 


8 . 5  Optical  Intensities  and  Electron  Densities 

The  purpose  for  showing  the  following  results  is  to  both  demonstrate 
a  working  daytime  ionospheric  chemistry-emission  code  and  to  provide  some  new 
information  on  the  use  of  optical  observations  for  deducing  electron  densities. 

We  shall  restrict  the  emission  results  to  the  LBH  1383  A  and  01  1356  A  features. 
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FLUX  I  /Cf1**2 


Figure  25  Comparison  of  calculated  and  measured  4tt 

integrated  Dhotoelectron  SDectra.  Both 
spectra  correspond  to  low  solar  activity, 
f q j  is  the  Jacchia  (1971)  n ( 0 )  scaling  factor 
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LOG|0  SOLAR  FLUX 
( photons/cm2  -  s  -  A) 


Figure  26  Solar  spectrum  for  low  solar  activity  used  in  the 
calculations.  Lines  have  been  given  a  reccangular 
distribution  with  a  5  A  width.  The  spectrum  be¬ 
tween  200  and  1000  A  is  from  Torr  et  al.  (1979)  who 
have  associated  with  it  a  10.7  cm  solar  flux  of  71. 

The  spectrum  below  200  A  come  from  Donnelly  and  Pope 

(1973).  ) 
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The  1383  A  feature  is  composed  of  the  (2,0)  and  (5,2)  bands.  The  LBH  band  at 
1354  A  contributes  only  a  sma1!  amount  to  the  total  1356  A  feature  and  for 
this  reason  will  not  be  considered  in  the  results  to  follow. 

Figure  28  shows  zenith  and  nadir  viewing  intensities  for  the  chosen 
features  under  three  sets  of  conditions.  Table  1  identifies  these  conditions 
which  refer  to  the  level  of  solar  activity  and  the  amount  of  atomic  oxygen  pre¬ 
sent.  The  case  numbers  in  the  table  are  used  to  label  the  curves  in  Figure  28. 

We  observe  a  change  in  intensity  of  2  in  going  from  low  to  high  solar  activity. 

As  expected,  reducing  the  0  density  by  a  factor  of  2  reduces  the  01  1356  A  intensity 
by  a  similar  factor.  The  LBH  1383  A  intensity  is  altered  only  slightly  because  of  a 
small  re-apportionment  of  photoelectron  energy  received  by  N~. 


Figure  29  shows  electron  densities  for  cases  1  and  2.  The  results  above 
200  km  include  the  effect  of  diffusion.  We  regard  these  results  as  preliminary 
since  they  are  the  first  obtained  by  code  CHEM  as  applied  to  the  daytime  problem 
and  since  the  diffusion  model  is  still  under  development.  As  preliminary  results, 
we  limit  them  to  cases  1  and  2  which  is  sufficient  for  the  demonstrating  variability 
in  the  electron  density  profile  and  how  it  relates  to  corresponding  variability  in 
optical  intensities. 
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4-1  (kR)  4rl  (kR) 


figure  28  Zenith  and  Nadir  viewing  01  1356  A  and  LBH  1383  A  dayglow 
intensities.  The  solar  zenith  angle  is  60°  and  Tm  =  1000°  K. 
See  Table  1  for  conditions  corresponding  to  the  three  cases 
considered. 
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Table  1 


Cases  for  which  results  have  been  obtained 
with  code  PEGFAC.  T  =  10^0°  and  the  solar 

oo 

zenith  angle  is  60° 


Case 

Jacchia  (1971) 
n(0)  scaling  factor 

F 

10.7 

1 

1 

71 

2 

1 

MB 

3 

.5 

u 

80 
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ABSTRACT 


In  this  paper  we  examine  the  relationship  among  certain 
prominent  auroral  FUV  emission  features,  the  incident  electron  spectrum, 
and  the  model  neutral  atmosphere.  Given  the  neutral  atmosphere  we  show 
that  for  simple  models  of  the  incident  electron  spectrum  (Maxwellian  and 
Gaussian  in  energy)  satellite  measurments  of  FUV  emission  features,  in 
principle,  determine  the  incident  electron  spectrum.  We  also  discuss  the 
relationship  between  the  incident  electron  spectrum  and  the  E-region 
plasma  density  profile  for  the  continuous  (diffuse)  aurora  and  for  a  stable 
arc. 
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INTRODUCT ION 


r' 


Our  not'ivation  for  studying  the  relationship  between  auroral 
UV  emissions  and  the  associated  electron  spectra  is  our  desire  to  util i 2e 
satellite  optical  emission  observations  to  infer  the  particle  spectra 
which  produce  them.  Under  certain  conditions  it  may  also  be  possible  to 
use  satellite  UV  measurements  to  infer  certain  E-region  ionospheric  proper¬ 
ties,  e.g.,  the  plasma  density  profile. 

This  is  the  first  of  two  papers  addressing  auroral  optical 
emissions  in  the  UV.  In  this  paper  we  examine  the  dependence  of  selected 
features  in  the  short  wavelength  region  of  the  Op  Schumann-Runqe  absorption 
continuum  on  the  incident  electron  spectrum.  This  region  lies  within  the 
FUV  region  which  runs  from  \  1050  A  to  1800  A.  As  will  be  seen,  the  varia¬ 
tion  cf  the  Op  absorption  with  wavelength  leads  to  some  noticeable  intensity- 
ratio  changes  as  the  hardness  of  the  incident  auroral  electron  spectrum 
changes.  The  ratio  we  refer  to  is  the  intensity  ratio  of  two  distinct 
emission  features  for  nadir  viewing  conditions. 

The  second  paper  now  in  preparation  will  address  features  in 
tne  midale  UV  such  as  01  2972  A  and  the  Vegard-Kapl an  bands.  As  the  wavelength 
increases,  the  question  arises  as  to  when  Rayleigh  scattering,  ground  and  cloud 
c’bedc,  end  scattering  from  aersols  become  important.  For  auroral  conditions 
the  arswer  is  at  wavelengths  longer  than  \  3000  A.  At  shorter  wavelengths  anc 
especially  for  the  region  to  be  addressed  in  this  paper,  pure  absorption  is 
strong  enough  to  prevent  the  radiation  from  penetrating  to  depths  where  scatter¬ 
ing  can  take  place. 

Electnons  and  protons  from  the  plasma  sheet  pitch-angle  scatter 
into  the  atmospheric  loss  cone  and  precipitate  to  produce  auroras.  In  the 
midnight  sector  both  electrons  and  protons  contribute  to  the  aurora  with  some¬ 
times  coincident  and  sometimes  separated  latitudina1  distributions  (Whalen 
and  Sharber,  1981).  The  proton  and  H-atom  contribution  to  the  auroral 
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ionization  rate  car  be  calculated  (  Jasperse  and  Basu,  198?)  and,  for  coin¬ 
cident  electron-proton  precipitation  and  typical  continuous  (diffuse,  auri>fc' 
conditions,  the  protons  and  H-atons  contribute  about  10  to  the  total  [-region 
ionization  rate  (Jasperse  and  Basu,  1982,  Figure  8).  In  this  paper  we  assume 
that  the  proton-H  atorr  contribution  to  the  FUV  emissions  is  also  about  10'  of 
the  electron  contribution  and  consider  only  the  effect  of  the  latter.  The  two 
types  cf  electron  produced  auroral  forms  we  shall  consider  are  the  continuous 
aurora  and  the  stable  arc.  They  are  discussed  in  detail  in  Section  3. 

There  have  been  many  attempts  at  deducing  the  incident  auroral 
spectrum  from  the  intensities  of  optical  emissions  in  the  visible  region  (see, 
e.g.  Rees  and  Luckey  (1974),  Vallance  Jones  (1975),  Arnoldy  and  Lewis  (1577), 
Shepherd  et  al.  (1980),  and  references  therein).  The  features  considered  have 
been  01  5577  A,  01  6300  A,  and  either  N£+  3914  A  or  4278  A,  all  of  which  are 
prominent  in  auroras  and  can  be  observed  from  the  ground.  For  most  analyses, 
the  intensity  variations  of  these  features  with  changing  electron  spectrum  com.e 
from  Rees  and  Luckey  who  characterized  their  spectra  by  Maxwellian  distributions. 

As  noted  by  Arnoldy  and  Lewis  (1977),  the  technique  has  been  fairly  successful 
in  estimating  the  mean  electron  energy  and  energy  content  of  the  incident  electron 
spectrum.  Rees  and  Luckey  do  acknowledge,  however,  a  potential  problem  due  to 
difficulties  over  the  years  in  explaining  the  5577  A  emission  which  is  chemistry 
dependent.  Examples  of  analyses  where  predictions  have  agreed  poorly  with  observa¬ 
tions  may  be  found  in  Rees  et  al.  (1977),  Sharp  et  al.  (1975),  and  Arnoldy  and 
Lewis  (1977).  For  this  reason  we  prefer  to  study  the  relationship  between  the 
electron  spectrum,  and  the  FUV  features  addressed  in  this  paper  which  have  the 
advantage  of  being  chemistry  independent. 


2.  BACKGROUND 


The  works  of  Strickland  et  al.  (1976),  Meier  et  al.  (1980), 

Meier  et  al.  (1981)  and  Conway  (1982)  form  the  basis  for  this  study.  The 
needed  electron  transport  model  comes  from  Strickland  et  al .  (1976).  Meier 
et  al .  (1980)  achieved,  for  the  first  time,  a  self-consistent  pictu-e  of 
the  FUV  h~  and  N3  dayglow  emissions  analyzing  the  rocket  data  of  Gentieu 
et  al .  (1979).  This  work  quantified  the  excitation  sources  and  processes 
affecting  the  transmission  of  the  photons,  namely  multiple  scattering  and 
cure  absorption.  Conway  (1982)  performed  a  line-by-line  synthesis  of  in¬ 
dividual  LBH  bands  examining  band  transmission.  He  derived  transmission 
functions  for  self  absorption  taking  into  account  the  varying  strengths  of 
the  individual  rotational  lines  and  found  the  transmission  to  be  larger 
than  previously  thought  (Meier  et  al.,  1981).  This  led  to  the  important 
conclusion  that  pure  absorption  dominates  self  absorption  for  all  of  the 

prominent  bands.  Meier  et  al .  (1981)  extended  their  dayglow  work  to  the 
analysis  of  auroral  rocket  data  of  Feldman  and  Gentieu  (1981)  and  again 
were  able  to  achieve  a  self-consistent  picture  for  the  above  features  as 
well  as  Nil  1085  A,  01  1356  A,  and  01  1304  A.  Regarding  LBH,  they  showed 
that  the  v'  populations  of  the  a  state  were  consistent  with  those  for 
direct  excitation  from  the  ground  state  .  They  were  also  able  to  infer  the  aor-- 
mate  nature  of  the  precipitating  electron  •''lux  and  the  amount  of  0  present  fro- 
the  available  LBH  anc  1356  A  data.  The  flux  determination  was  made  possible  b.> 
the  differing  amounts  of  0^  pure  absorption  from  one  LBH  band  to  another  which 
is  the  effect  to  be  discussed  in  some  detail  in  this  paper. 


We  are  considering  the  following  features:  IN  (0,0)  at 

3914  A,  03  1356  A,  and  the  N?  LBH  bands  at  1273  A  (6,0),  1325  A  (4,0), 

1354  A  (3,0),  and  1384  A  (2,0  and  5,2).  Spectroscopi cally ,  3914  A 

stands  apart  from  the  rest  of  the  features  but  is  included  since  its  emis¬ 
sion  is  strong,  it  is  measured  in  most  auroral  particle  and  optical  experi 
ments,  energy  input  relates  directly  to  its  intensity,  and  finally  it  is 
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well  unaerstooc.  Tr<c  rem«nring  features  hove  beer,  selectee  i  1 '  a-  1)  they 
possess  no  albedo  ( r  or  below  the  ionosphere,  2)  some  show  strong  variations 
with  electron  spectral  hardness,  3)  all  are  prominent  in  auroral  FUV  spectra, 

4)  blending  is  not  a  problem  since  they  do  not  blend  with  one  another  (with  the 
exception  of  01  1156  4  and  LBH  1354  A  as  noted  below)  nor  with  other  nearby 
features  such  as  01  1304  A  for  resolutions  better  than,  say,  10  A  (see,  e.g., 
Feldman  and  Gentieu,  1981),  and  finally,  5)  these  features,  like  3914  A, 
are  reasonably  well  understood  (Meier  et  al.,  1980,  Meier  et  a!.,  1961).  LBH 
3354  A  and  01  1356  A  are  effectively  blended  and  for  this  reason  will  be  discusseo 
at  times  in  this  paper  as  one  feature.  We  do  not  wish  to  eliminate  them  from  con¬ 
sideration  because  of  blending  since  01  1356  A  has  the  potential  as  a  useful 
signature  on  the  amount  of  atomic  oxygen  (0)  present.  Mention  should  be  made 
of  01  1304  A  since  it  Is  one  of  the  most  prominent  features  in  the  FUV  spectral 
region.  We  choose  not  to  include  this  feature  in  the  above  list  for  two  reasons. 
First,  there  is  a  severe  computational  problem  in  accurately  determining  its 
emission  properties  due  to  extensive  multiple  scattering  (Strickland  and  Donahue, 
1970;  Meier  and  Lee,  1982).  Second,  its  nadir  intensity  is  not  necessarily  an 
accurate  signature  of  electron  impact  excitation  along  the  given  line  of  sight 
(even  correcting  for  the  obvious  multiple  scattering  effects)  due  to  the  non¬ 
local  nature  of  the  given  emission.  If,  for  example,  a  nadir  observation  were 
made  adjacent  to  some  region  of  brighter  emission,  some  of  the  observed  signal 
would  be  due  to  radiation  originating  from  this  brighter  region  and  would  con¬ 
sequently  overestimate  the  initial  electron  impact  excitation  in  the  observed 
rec  i  or. . 

Zipf  and  his  colleagues  (McLaughlin  et  al.,  19821  have  recently  reported 
that  various  LBH  bands  are  blended  with  atomic  nitrogen  lines  and  other  bands. 

In  addition,  they  suggest  that  optical  pumping  of  the  a  g  state  from  the 
c^'*iu+  state  is  a  significant  source  in  the  upper  atmosphere.  These  complications 
will  not  be  addressed  in  this  paper.  First,  such  complications,  from  a  compu¬ 
tational  point  of  view,  will  only  affect  the  applied  emission  cross  section  for 
each  feature  so  affected.  This  could  lead  to  a  change  in  the  strengths  of  what 
we  are  referring  to  as  our  LBH  features  but  not  any  noteable  changes  in  their 
relative  behavior  from  one  auroral  situation  to  another.  This  last  statement 


applies  because  the  same  impact  species  (Np )  1S  responsible  for  the  mulnplc- 
components  suggested  by  McLaughlin  et  al.  Second,  our  interest  is  in  optical 
properties  on  the  resolution  scale  appropriate  to  today's  rocket  and  satellite 
borne  UV  spectrometers  which  do  not  resolve  blended  features  as  discussed  above. 

We  shall  continue  to  refer  to  the  features  being  considered  as  was  done  in  the 
previous  paragraph  keeping  in  mind  that  some  scaling  may  be  called  for  in  the  future. 

We  noted  above  that  there  should  be  little  relative  intensity  variation 
among  those  features  being  referred  to  as  LBH  bands  due  to  either  blending  or 
cascading  as  auroral  conditions  change.  This  statement  is  based  on  the  assump¬ 
tion  that  in  going  from  one  auroral  situation  to  another,  the  electron  spectrum 
has  not  shifted  from  its  expected  occurrence  in  the  keV  region  to  the  sub-keV 
region.  As  long  as  we  exclude  that  situation,  there  is  little  relative  variation 
in  excitation  efficiencies  between,  say,  the  a  \l  and  c'  states  as  the  electron 
spectrum  is  changed.  One  does  arrive  at  a  point,  however,  as  the  spectrum  softens, 
where  this  is  no  longer  true  due  to  the  differences  in  respective  impact  cross 
section  shapes.  For  very  soft  auroral  spectra,  the  more  forbidden  nature  of  the 
a  cross  section  leads  to  relatively  more  a  \l  excitation  compared  to,  say,  c1 
excitation.  Our  transport  calculations  show  that  this  effect  is  not  important 
as  long  as  the  mean  energy  of  the  incident  electrons  is  greater  than  approxi¬ 
mately  one  kilovolt. 

There  is  one  obvious  disadvantage  to  using  the  chosen  set  of 
features  for  probing  the  auroral  ionosphere  in  contrast  tc  the  set 
3514  A,  01  5577  A,  and  01  6300  A  which  we  discussed  earlie1'.  That  dis¬ 
advantage  concerns  the  demands  placed  on  the  measurements .  A  spectrometer 
is  reguired  for  which  one  must  be  concerned  with  resolution,  scanning  rates, 
and  signals  weaker  than  from  such  prominent  features  as  those  just  cited. 

This  point  is  further  emphasized  by  the  present  paucity  of  good  spectral  data 
from  satellite  experiments  for  our  chosen  set  of  features.  Nevertheless, 
demands  placed  on  instrumentation  are  becoming  more  modest  with  present 
advancements  in  state-of-the-art  detection  systems. 
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3.  SOME  FUV  EMISSION  OBSERVATIONS 


An  auroral  optical  spectrum  including  the  FUV  features  o f 
interest  to  this  work  is  given  in  Figure  1.  This  spectrum  was  obtained 
in  a  rocket  experiment  with  the  spectrometer  looking  down  from  above  the 
emitting  layer  (P.D.  Feldman,  private  communication,  1981).  The  02  Pure 
absorption  cross  section  is  included  as  well  as  an  altitude  scale  showina 
where  unit  optical  depth  occurs  due  to  02  absorption  based  on  the  adjacent 
cross  section  scale.  In  Section  6,  we  will  show  that  absorption  strongly 
affects  the  LBH  band  at  1384  A  but  has  little  effect  on  the  band  at  1373  A 
over  the  range  of  incident  electron  spectra  considered. 

There  are  two  basic  types  of  E-region  auroras;  continuous  (or 
diffuse)  and  discrete.  The  two  are  different  not  only  in  spatial  character¬ 
istics  but  also  in  temporal  and  precipitating  particle  characteristics.  The 
two  types  are  illustrated  in  Figure  2  which  shows  a  latitudinal  graph  of  a 
band  of  FUV  emissions  detected  by  a  nadir  viewing  photometer  on  satellite  S3-4 
as  it  crossed  the  auroral  region  near  midnight  (Huffman  et  a  1 . ,  1981).  The  band 
pass  (full  width  at  half  maximum)  of  the  filter  used  was  165  A  with  its  trans¬ 
mission  maximum  at  1550  A.  The  broad  continuum  which  peaks  near  -  66. 5C  C.G. 
latitude  shows  the  continuous  aurora;  the  peaks  near  -  69c  and  -  70. 5r  identify 
discrete  auroras. 

The  latitudinal  distribution  of  the  continuous  auro-a  has  an 
FWHM  of  about  3°  and  so  is  consistent  with  other  observations  which  find  the 
distribution  of  precipitated  particle  energy  flux  to  be  near  Gaussian  in  latitude 
with  approximately  the  same  scale  parameter  (Whalen,  1981).  These  observations  also 
find  the  latitudinal  distribution  to  be  uniform  throughout  approximately  12  hours 
of  magnetic  local  time.  The  distribution  in  Figure  2  departs  somewhat  from  a 
Gaussian  because  the  FUV  emission  is  dependent  not  only  on  energy  flux,  but  also 
on  particle  characteristic  energy  which  is  known  to  vary  with  latitudes  (Sharber, 
1981). 
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Although  a  very  extensive  and  stable  auroral  feature,  the 
continuous  aurora  has  frequently  not  been  detected  because  of  its  relatively 
low  brightness.  Even  when  detected,  the  morphology  of  this  aurora  has  beer 
difficult  to  describe  because  its  distribution  broadens  in  latitudinal  width 
with  decreasing  intensity.  As  a  result,  its  apparent  width  is  a  function  of 
the  threshold  of  the  detector  which  observes  it.  Thus,  the  distribution  of 
continuous  aurora  in  Figure  2  would  appear  to  be  about  2C  wide  to  a  detector 
with  a  threshold  corresponding  to  100  R,  but  approximately  8C  wide  to  a  detector 
with  a  threshold  of  6  R.  Similar  difficulties  exist  in  the  specification  of 
the  latitude  of  the  equatorward  boundary  of  this  aurora  because  of  its  sensitive 
dependence  on  detector  threshold. 

The  continuous  aurora  is  stable  on  a  time  scale  of  at  least 
minutes  (Whalen,  1981)  and  is  produced  by  precipitating  keV  electrons  with  a 
minor  contribution  from  protons  (Lui  et  al.,  1977;  Sharber,  1981).  The  energy 
spectrum  of  the  electrons  obtained  by  averaging  over  the  loss  cone  is  usually 
Maxwellian  (Sharber,  1981). 

Discrete  auroras,  on  the  other  hand,  are  extremely  variable 
in  time,  latitude,  and  longitude.  The  normal  behavior  of  the  discrete  form 
is  to  change  on  a  time  scale  that  is  short  compared  to  an  E-  region  chemicai 
time  constant  ("■  1  minute).  Discrete  auroral  forms  that  vary  on  a  tine  scale 
longer  than  a  chemical  time  constant  will  be  regarded  as  stable  arcs.  The 
prominent  auroral  arc  in  Figure  2  was  imaged  by  a  DMSP  satellite  at  aporo>i- 
mately  the  same  location  at  which  the  S3-4  satellite  observed  it  (Huffman  et  al . , 
1S81).  Since  the  DMSP  measurement  preceded  that  of  the  S3-4  satellite  by 
10  minutes,  the  arc  was  apparently  stable  in  position.  However,  there  is 
no  information  which  could  insure  that  the  arc  was  also  stable  in  brightness. 
Discrete  auroras  are  produced  by  precipitating  electrons  which  have  energy 
spectra  that  are  often  modelled  by  Gaussian  functions  of  energy  (e.g. 

Fontheim  et  al . ,  (1982) 
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The  source  of  particles  that  produce  both  types  of  auroral 
forms  is  apparently  the  plasma  sheet;  those  electrons  producing  the  continuous 
aurora  usually  have  Maxwellian  spectra  characteristic  of  the  plasma  sheet 
(Rearwin  and  Hones,  1974)  and  those  electrons  producing  discrete  auroras  have 
spectra  consistent  with  additional  acceleration  produced  by  electric  field 
structures  located  in  the  magnetosphere  (Evans,  1974).  In  this  regard  the 
maximum  in  the  near  Gaussian  latitudinal  distribution  for  the  continuous 
aurora  is  an  important  magnetospheric  feature  since  it  represents  the  normal 
equatorward  limit  of  discrete  auroras  (Lui  et  al . ,  1977).  In  Figure  2  we  note  the 
absence  of  discrete  auroras  in  the  equatorward  half  of  the  latitudinal  distribution. 
Furthermore,  because  of  the  uniformity  of  the  continuous  aurora  in  local  time, 
the  value  of  the  energy  flux  at  the  maximum  in  the  latitudinal  distribution  is  an 
indicator  of  auroral  activity  on  a  large  scale  (Whalen,  1981).  A  further  fundamental 
significance  of  the  maximum  in  the  latitudinal  distribution  shown  in  Figure  2  is 
its  location  at  or  near  the  interface  between  positive  and  negative  field  aligned 
currents  (Robinson  et  al.,  1982). 


95 


4.  APPROACH 


For  the  chosen  set  of  emission  features  we  wish  to  know  their 
relative  variations  in  intensity  as  functions  of  the  incident  electron  spectrum 
and  the  model  atmosphere.  To  study  this  problem  we  solve  the  Boltzmann  electron 
transport  equation  for  several  incident  differential  electron  number  fluxes 
characterized  by  Maxwellian  and  Gaussian  distributions.  The  equation  is 


P do ( z , E ,  u ) 
dz 


k(z,E)  [  -  c(z,E,y) 


r2(e'  ,E,e)o(z,E'  ,U'  )d.r  dE 1  ] 


+  V2)  fr  £lp(z,e)4>(2.e,u)] 


with  terms  defined  as  follows: 

z,E,u  altitude,  energy,  and  direction  cosine 

?  scattering  angTe 

c  differential  electron  number  flux  (electron  spectrum)  in 

2 

e/cm  -s-eV-sr 

k  total  inverse  mean  free  path  (1MFP)  in  cm" ‘ 

r ,  fractional  IMFP  for  species 

Rj.  scattering  and  energy  loss  redistribution  function  for 

species  in  eV^sr"1 

np  plasma  density  in  cm'^ 

2 

Lp  plasma  loss  function  in  eV-cm  . 
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The  function  k  is  given  by 


k(z,E)  = 


n£(2)o£(E) 


[cm'1] 


(?) 


where  n.  is  the  neutral  density  (Ng.Og,  or  0)  and  z;  is  the  total  electron 
impact  cross  section  (elastic  plus  inelastic).  The  function  r.  is 

r,(z,E)  =  [k(z,E)r1n.(z)c;(E)  (3) 

The  function  R  contains  several  terms  describing  elastic  scattering,  excitation, 
and  ionization.  More  details  regarding  R£  and  the  solution  of  Equation  (1)  mav  be 
found  in  Strickland  et  al.  (1976). 

Equation  (1)  yields  a  solution  from  primary  electron  energies  in 
the  keV  range  to  secondary  electron  energies  in  the  low  eV  range.  Thus,  the 
volume  emission  rates  of  interest  to  us  here,  lacking  chemically  produced 
components,  may  be  obtained  from  this  solution  using 


P-,lz) 

K  I 


Cjj.^  (E)c(z  ,E,v)2-dudE 


r  -3  -1, 

[cm  -s  ] 


(4) 


where  wfci  is  the  threshold  for  the  ith  excitation  process  of  the  Ith  species. 
The  corresponding  intensity  in  the  absence  of  enhanced  emission  caused  by 
multiple  scattering  is 
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The  function  k  is  given  by 


k(*,E)  =  Yi  [cm'1] 

K 


where  nf  is  the  rt'h  neutral  density  (N2 .02 »  or  0)  and  c.  is  the  total  electron 
impact  cross  section  (elastic  plus  inelastic).  The  function  r.  is 


r .  (2  ,E)  =  [k(z,E)]_1n.(z)c,(E) 

x.  >.  >■ 

The  function  R  contains  several  terms  describing  elastic  scattering,  excitation, 
and  ionization.  More  details  regarding  R£  and  the  solution  of  Equation  (1)  may  be 
found  in  Strickland  et  al.  (1976). 

Equation  (1)  yields  a  solution  from  primary  electron  energies  ir 
the  keV  range  to  secondary  electron  energies  in  the  low  eV  range.  Thus,  the 
volume  emission  rates  of  interest  to  us  here,  lacking  chemically  produced 
components,  may  be  obtained  from  this  solution  using 


P-Az) 

K  ^ 


c. .  (E)tU  ,E,p)2rdudE 

A.  1 


»•'] 


where  .  is  the  threshold  for  the  i*11  excitation  process  of  the  I**1  species. 
The  corresponding  intensity  in  the  absence  of  enhanced  emission  caused  by 
multiple  scattering  is 
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':VZ-V  ■  {“- 


-r'f  pri  <z'  «(>[- ;  (tJ,((z'  )-t;  (  U) 0  jT(z'  .  0 ) cr  ■ .  ( 


n  i  n 


[ph/cn-  -  S-sr] 


where  ;-0  is  the  cosine  of  the  look  direction  (uQ  =  1  for  true  nadir  direction),  t  is 
the  pure  absorption  optical  depth,  and  T  is  the  self  absorption  transmission  func¬ 
tion  (see,  e.g.  Strickland,  1979  for  lines  and  Conway,  1981  for  bands).  rhe  factor 


Wj .  is  a  branching  ratio  and  allows  for  the  possibility  of  more  than  one  emission 
feature  arisinq  from  the  £ith  excited  state.  For  example,  the  index  j  for 
molecular  bands  represents  the  double  index  v'  v"  designating  the  upper  and  lower 
vibrational  levels  of  the  transition  producing  Jfie  band  emission  pesults  to 

follow  will  be  in  the  form  of  4ttI  with  units  of  kilo-Rayleighs.  shall 

continue  to  use  the  term  intensity  for  this  quantity  although  mo  nroperly 
it  is  referred  to  as  an  apparent  column  emission  rate.  The  dept  ' n  given 
for  and  has  the  form 


t(z)  = 


nn  (z')dz' 
u2 


figure  3  shows  t  versus  z  for  the  LBH  bands  under  investigation.  We  note  that 
unit  optica1  depth  occurs  at  less  than  90  km  for  the  1273  A  feature  in  contrast  r 
120  km  for  the  1384  A  feature.  In  turn,  as  stated  earlier,  the  1273  A  intens't.v 
is  not  affected  by  pure  absorption  over  the  range  of  incident  electron  spectra 
being  considered  whereas  the  1384  A  intensity  is  strongly  affected. 

A  modified  ve'rsion  of  Equation  (5)  should  be  used  for  01  1356  A 
due  to  multiple  scattering  effects  (Anderson  et  al.,  1980).  Its  form  is 


I(z,u0)  =  (4m)'1 y* 


S(z’)  [f  T5(z',z,u0)  +  |  T3(z',z,u0)]dz7.( 


min 
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w here  S  is  t r. t  tota'  volume  emission  rate  ir,  both  lines  of  this  doublet  tu't 
(  1355.5  A  and  1356.5  A;  .  The  transmission  fraction  refers  to  1  355.5  A  v»r  :  1  e 

its  subscript  refers  to  the  statistical  weight  of  the  ^  ground  state  le,e'. 

The  function  S  is  obtained  from  the  integral  equation 

/2max  r  3 

S ( 2 1 )H(z ,z‘ )dz‘  [ph/cm  -s]  (£ 

zmin 

where  SQ  is  the  initial  excitation  rate  [given  here  as  P  .  in  Equation  (4)]  and 
H  is  a  modified  Holstein  H  function  (Holstein,  1947)  taking  into  account  the 
doublet  structure  of  the  1356  A  feature.  A  detailed  discussion  of  Eouaticn  {£, 
as  it  applies  to  another  multiplet,  that  of  01  1304  A,  is  given  by  Strickland 
and  Donahue  (1970).  The  differences  arising  in  4rl  by  using  Equation  (7)  in 
place  of  Equation  (5)  will  be  discussed  later. 

We  have  chosen  to  describe  some  of  our  incident  electror  spectra 
with  Maxwellian  distributions  since  they  are  characteristic  of  the  precipitation 
producing  the  continuous  aurora  as  determined  by  averaging  spectra  over  the 
loss  cone  (Sharber,  1981  and  Whalen  and  Sharber,  1981).  In  addition,  Maxwellians 
are  convenient  to  use  and  have  been  frequently  used  before  (see,  e.g.,  Roble  and 
Rees,  1977  and  Rees  and  Luckey,  1974).  Let  dQ(E)  be  the  incident  differential 
electron  number  flux  which  we  assume  to  be  isotropic.  For  a  Maxwellian,  it  is 
given  by 

♦0(E)  =  (~^4)  E  exp  [-  E/E*]  [e/cm2-s-eV-sr] 

'2tEm 

where  E  is  the  electron  energy,  E^  is  the  characteristic  energy  and  is  the 
energy  flux.  Under  certain  conditions  the  incident  spectrum  includes  a  high-enerov 
tail  (Vondrak,  1981).  A  high-energy  tail  can  be  included  by  adding  a  second  term 
to  Equation  (9)  and  would  introduce  additional  parameters.  Since  the  high-energy 
tail  is  not  always  present  and  since  it  produces  mainly  D-region  effects  we 
neglect  its  effect. 
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Maxwellian  distr  ibutions  dc  not  characterize  au'-pral  a>  c  e^ec’-'O"  spectra 
very  satisfactorily  since  these  spectra  generally  have  a  s^arof  pea*  .'tee. 

Bo;, d  ,  1575  and  Lui  et  a1..,  1977  ).  For  this  reason  we  have  included  so-c  ru'-rc* 

Gaussian  distributions  ir  our  study.  The  expression  for  :  using  a  Ge-ssien 
distribution  is 


'o 


exp 


[e/crr  -S*eV-Sr] 


(10) 


where  the  quantities  introduced  have  similar  meanings  as  those  introduced  ir, 
Eauation  (9).  In  this  work,  the  Gaussian  scale  parameter,  W,  has  been  c’ven 
values  equal  to  0.25  Eg.  In  the  following  sections,  when  general  properties  of 
the  results  are  discussed,  we  shall  use  E„  and  Q„  to  refer  to  either  EM  or  Er 

0  0  no 

and  to  Qm  or  Qg. 
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5.  MODEL  ATMOSPHERE  AND  ATOMIC  AND  MOLECULAR  PARAMETERS 


The  basic  model  atmosphere  'sed  in  this  work  comes  from  Jacchia 
(1977).  This  is  given  in  Table  1.  An  indication  of  density  variability  in 
the  auroral  ionosphere  is  given  in  Figure  4  where  the  Jacchia  model  is  com¬ 
pared  with  one  reported  by  Rees  et  al.  (1977)  and  Sharp  et  al .  (1979).  The 
comparison  is  orovided  through  ratios  of  the  respective  model  densities. 

Atomic  oxygen  is  seen  to  be  more  than  a  factor  of  three  smaller  as  reported 
by  Rees  et  al .  and  Sharp  et  al .  Guided  by  such  differences,  some  o*  our  re¬ 
sults  have  teen  obtained  with  a  scaled  down  Jacchia  0  density.  Specifically, 
intensities  will  be  shown  for  the  Jacchia  0  density  reduced  by  three.  There 
is  also  0^  density  variability  shown  in  Figure  4  although  this  is  modest  below  12 
where  pure  absorption  exhibits  most  of  its  effect.  Results  have  beer  obtained 
with  two  scalings  of  the  Jacchia  density.  The  scaling  factors  are  l.E  and 
.667  which  allow  for  a  range  of  variables  that  we  assume  is  more  thar  adecate 
to  reflect  actual  changes  below  120  km. 

A  full  set  of  electron  impact  ionization  and  excitation  cress  sections 
for  N0,  0^  and  0  is  needed  to  perform  the  transport  calculations.  The  basic  set 
is  given  by  Strickland  et  al .  (1976)  and  Oran  and  Strickland  (1978)  with  so~e 

modifications  based  on  recent  experimental  data.  The  needed  rates  c.  reqL're 

1  5  1  ‘  + ,  2  ♦ 

cross  sections  for  exciting  N^a  *"  '  and  71  (  S)  and  icczirg  N,  tc  N„  IE 

C  C  c  £ 

~at  l e  Z  gives  peek  cess  sectic  values  •  arc  references  tc  tneir 

full  energy  dependences .  "nese  are  t he  sane  coss  sections  used  in  the  aurora' 

analysis  of  Meier  et  al.  (1981).  Table  3  gives  the  percent  of  N?(a  )  and 
-f  2  +  ^9 

N^  (B  Du  )  production  going  to  the  various  bands  under  consideration.  The 

table  also  contains  the  0^  Pure  absorption  cross  sections  for  all  of  the 

features.  Again,  these  values  are  the  same  as  those  appearing  in  Meier  et  al. 

(1981).  The  emission  factor  in  Table  3  for  3914  A  cones  from  Vallance  Jones 

(1974).  The  factors  for  LBH  were  derived  by  R.R.  Conway  and  are  similar  to 

those  given  by  Vallance  Jones  for  the  specific  bands  of  interest. 
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Ar.  assumption  we  make  in  this  work  is  that  the  v'  populations  of  the 
a  1  .  state  do  not  vary  with  changing  auroral  conditions.  In  turn,  this  also  "tsrs 
that  the  emission  factors  do  not  change.  This  assumption  should  be  valic  as 

lorg  as  additional  excitation  mechanisms  beyond  those  discussed  in  Secfinr  2 
are  not  operating.  Huffman  et  al.  (19B0)  have  reported  satellite  photometric 
and  spectrometer  observations  of  nighttime  LBH  emissions  which  extend  into  the 
auroral  region.  The  nighttime  band  distribution  is  shifted  to  longer  wavelengths 
compared  to  that  for  auroral  electron  impact  excitation.  Such  a  nighttime  component 
of  a  low  level  auroral  spectrum  may  change  the  assumed  v'  distribution  since  total 
nighttime  LBH  intensities  in  excess  of  1  kR  have  been  observed.  Caution  should, 
therefore,  be  exercised  in  using  results  to  follow  for  weak  auroras. 
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6.  CALCULATED  VOLUME  PRODUCTION  RATES  AND  INTENSITIES 


6.1  Volume  Production  Rates  as  Functions  of  Eq 


Results  have  been  obtained  for  several  values  of  E.,  and  E,~. 

n  b 

These  are  EM  =  .5,  1.0,  2.5,  and  5.0  keV  and  Eg  =  1.0,  2.0,  3.3,  5.0,  7.5 
and  10.0  keV.  The  energy  parameters  and  Qg  have  been  selected  to  give 
1  erg/cm2-s  for  the  incident  energy.  Figure  5  gives  the  N?  a  production 
rates  for  the  incident  spectra  given  by  the  above  characteristic  energies. 
The  top  panel  contains  the  Maxwellian  results  while  those  for  the  Gaussian 
spectra  appear  in  the  lower  panel. 


As  expected,  the  Gaussian  results  are  more  sharply  peaked 
and  one  cannot  find  an  Eg  value  which  precisely  reproduces  the  re¬ 
sults  for  some  given  EM  value.  The  relationship  Eg  =  3EM  is  approximately 
correct,  however,  for  achieving  the  same  altitudes  of  maximum  excitation 
and  ionization. 

Figure  6  shews  production  rates  for  0  (5S).  For  the  larger 
characteristic  energies,  the  area  beneath  a  given  profile  (column  production 
rate)  decreases  with  increasing  characteristic  energy  in  contrast  with  the 
behavior  in  Figure  5.  This  is  caused  by  the  decrease  in  the  0  to  N^  density 
ratio  with  decreasing  altitude  and  has  the  desirable  effects  for  producing 
stronger  intensity  variations  with  electron  spectral  hardness  than  would 
occur  by  pure  absorption  effects  alone. 


6-2  Intensities  as  Functions  of  Eq  and  Qq 

For  our  standard  model  atmosphere  (Table  1),  intensities 
appear  in  Figures  7  and  8  versus  characteristic  energy  for  Maxwellian  and 
Gaussian  incident  electron  spectra.  The  dominant  effect  producing  the 
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differing  variations  with  £q  among  these  intensities  is  pure  absorption. 

At  small  values  of  Eq,  pure  absorption  has  litlle  effect  on  any  of  the  fee* 
tures  since  emission  occurs  where  t^  is  small.  The  fall  off  observed  in 
the  emissions  as  Eq  approaches  ze£o  results  from  high  altitude  energy 
deposition  where  there  is  relatively  less  compared  to  0  than  at  lower 
altitudes.  Consequently,  less  of  the  available  enerov  is  beinc  absorbed  by 
and  more  by  0  which  results  in  weaker  Ng  emissions.  At  lc''ger  values  o{ 

Eq,  most  emission  occurs  where  strong  contrasts  exist  in  pure  absorption 
from  one  feature  to  another.  We  see  that  this  leads  to  relative  variations 
as  large  as  a  factor  of  about  four  between  1273  A  and  1384  A  for  the  given  Eq 
range.  Variations  of  this  strength  make  such  features  attractive  as  monitors 
of  the  hardness  of  the  incident  electron  spectrum.  The  weakness  of  their 
intensities  relative  to  3914  A  may  make  it  difficult  to  obtain  good 
signal-to-noise  ratios  in  actual  observations  unless  the  incident  spectrum 
contains  several  ergs/cm  -s.  It  may  be  more  attractive  to  monitor  the 
intensity  at  3914  A  in  place  of  1273  A  due  to  its  brightness,  its  similar 
variation  with  Eq,  and  the  ability  to  make  the  measurement  with  a  photometer 
although  account  must  be  taken  of  albedo  and  backscatter  unlike  in  the  FUV 
region . 


The  strongest  variations  with  Eq  occur  for  01  1356  A  end  result 
from  the  combined  effects  of  0^  pure  absorption  and  the  relative  variation 
of  0  to  N'2  as  a  function  of  altitude.  For  larger  values  of  Eq  .  'ess  of  the 
available  incident  electron  energy  goes  to  0  due  tc  lack  o-r  th-s  constituent 
where  the  energy  deposition  i<  occuring.  01  1356  A  is  cleanly  a  useful  feature 
for  estimating  electron  spectral  hardness  if  one  knows  the  0  density.  We  have 
seen  evidence  of  enough  variability  in  this  density  that  such  usefullness  can 
be  expected  to  be  rather  limited  without  some  measure  of  this  quantity.  Sub¬ 
section  6,3  will  address  the  sensitivity  of  the  1356  A  intensity  to  scalings 
of  the  0  density. 

We  noted  earlier  that  there  is  some  effect  on  the  01  1356  A 
intensity  due  to  multiple  scattering.  The  extent  of  this  effect  has  been 
determined  for  selected  Eq  values  and  two  n(0)  profiles  (D.E.  Anderson,  Jr., 
private  communication,  1982).  One  profile  is  that  in  Table  3,  while  the 
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other  is  obtained  by  scaling  the  first  by  .33.  The  difference  between  the 
single  scattering  (Equation  5)  and  multiple  scattering  (Equation  7)  intensi¬ 
ties  is  about  15*  at  small  Eg  values  for  n(0)  given  in  Table  3.  It  decreases 
with  increasing  Eg  because  of  the  increasing  importance  of  0^  pure  absorption 
which  inhibits  multiple  scattering.  For  the  reduced  n(0)  profile,  the  maximum 
difference  is  only  about  5%.  In  either  case,  multiple  scattering  effects  will 
not  significantly  impact  upon  our  long  range  goals  in  this  work.  Tney  do  be¬ 
come  more  important  when  viewing  along  greater  1 ines-of-si ght  and  should  not 
be  ignored,  e.g.,  when  considering  limb  viewing  conditions  (Strickland  and 
Anderson,  1983). 

Fiqure  9  shows  the  intensities  as  functions  of  Qg  for  the  Max¬ 
wellian  characteri stic  energy  E^  =  2.5  keV.  This  value  was  selected 
because  it  lies  in  the  middle  of  the  applied  range.  Since  the  applied 
Boltzmann  equation  is  linear,  its  solution  as  well  as  the  intensities 
based  upon  it  vary  linearly  with  Qg.  In  spite  of  this  simple  behavior, 

Figure  9  is  included  since  it  is  useful  for  obtaining  the  various  emission 
strengths  at  a  glance  for  a  range  of  the  precipitating  electrons'  energy 
content. 

6.3  Intensities  as  a  Function  of  the  Altitude  for  Maximum  Energy  Deposition 

It  is  instructive  to  display  the  intensities  in  Figures  7  and  6  as 
functions  of  the  altitude  of  peak  energy  deposition,  2  ,  which  bears  a  one-to-one 
correspondence  with  Eg.  By  doing  so,  we  may  combine  the  intensities  for  Max¬ 
wellian  and  Gaussian  electron  spectra  on  the  same  plot  and  conveniently  note 
similarities  and  differences  in  the  intensity  behavior  for  these  two  representa¬ 
tions. 


Figure  10  shows  the  intensities  of  Figures  7  and  8  versus  the  given 
altitude  parameter.  For  convenience,  let  us  refer  to  the  intensities  for  the 
two  source  representati ves  as  1^  and  Ig.  In  general,  lg  decreases  more  slowly 
with  decreasing  2  .  To  understand  why,  we  must  consider  the  effect  of  pure 
absorption  for  volume  emission  profiles  peaking  at  the  same  altitude  but  with 
different  distributions  about  this  peak.  It  is  clear  that  a  smaller  amount 
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of  the  emission  below  the  peak  is  purely  absorbed  for  the  Gaussian  case  since 
the  emission  does  not  extend  as  deeply  into  the  atmosphere.  The  opposite  effect 
occurs  above  the  peak  where  relatively  more  of  the  emission  for  the  Maxwellian 
case  extends  into  the  non-absorbi ng  region.  Pure  absorption  under  these  circum¬ 
stances  may  lead  to  either  greater  or  less  intensity  for  one  type  of  intensity 
versus  the  other.  Greater  intensities  have  been  seen  to  occur  in  Figure  ID  for 
the  Gaussian  case  over  the  range  of  Eg  considered.  Less  intensity  will  finally 
occur,  however,  with  increasing  Eg  since  pure  absorption  becomes  sc  strong  that 
only  the  upper  part  of  the  volume  emission  profile  contributes  to  the  intensity. 
Since  this  profile  has  relatively  less  extension  to  high  altitudes  compared  to 
the  Maxwellian  case,  relatively  less  radiation  can  escape.  The  altitude,  z 
at  which  this  occurs  is,  of  course,  a  function  of  the  feature  or  wavelengtr. 

For  all  features  shown,  this  altitude  occurs  below  104  km. 


Our  primary  motivation  for  generating  Figure  10  was  tc  determine 
how  well  Zp  can  be  specified  from  intensity  ratios  regardless  of  electron  source 
characteristics.  Clearly,  these  ratios  are  not  useful  for  soft  sources  giving 
Zp  >'  115  km.  The  situation  is  much  improved  for  harder  sources  leading  to  the 
exhibited  variations.  A  unit  ratio  of  1 1384  t0  ^1273*  ^or  examp"le’  9 1 v e s 
a  value  close  to  105  km  for  either  representation.  Even  better  determination 
of  z  is  possible  using  4-Ij,^  if  one  knows  the‘0  density  profile. 


6.4  Intensities  as  a  Function  o'  the  Model  Neutral  Atmospnere 

We  introduced  our  applied  model  atmosphere  in  Section  5  and  noted 
the  kinds  of  density  variations  that  may  occur  under  varying  auroral  conditions. 
Uncertainties  in  the  0^  and  0  densities  can  have  a  strong  impact  on  the  use- 
ability  of  nadir  intensities  for  the  features  being  considered.  For  this  reason, 
we  have  determined  the  sensitivity  of  these  intensities  to  what  we  regard  as 
reasonable  ranges  of  variability  in  the  above  densities.  Figure  11  gives  the 
1356  A  intensity  versus  EM  for  the  Jacchia  0  density  and  one  scaled  down  by  a 
factor  of  3.  The  LBH  contribution  is  unaffected  by  this  variation.  The  01  1356  A 
component,  on  the  other  hand,  varies  almost  linearly  with  scalings  of  n(0).  This 
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behavior  clearly  limits  the  usefulness  of  the  1356  A  intensity  as  an  indicator 
of  electron  spectral  hardness  unless  n(0)  is  known.  We  may  turn  the  argument 
around,  however,  and  note  that  if  the  hardness  is  known  by  other  means,  then 
this  feature  will  be  useful  for  determining  the  amount  of  0  present. 

We  now  consider  the  sensitivity  of  selected  intensities  to 
n^)  variations.  The  results  appear  in  Figure  12  for  n ( 0^ )  ranging  from  .67 
to  1.5  times  the  Jacchia  values.  The  given  variations  clearly  complicate 
matters  if  nfO^)  cannot  be  better  specified  than,  say,  50%  in  actual  calcula¬ 
tions.  Specifically,  we  are  addressing  uncertanities  below  n,  120  km.  At  these 
lower  altitudes,  the  model  comparisons  in  Figure  4  lead  to  much  smaller  dif¬ 
ferences  which  if  indicative  of  the  actual  degree  of  n ( 0^ )  variability  in  the 
auroral  ionosphere  suggest  that  such  variability  will  not  play  an  injurious 
role  in  the  use  of  intensities  as  proposed  here. 
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7.  CALCULATED  IONIZATION  RATES  AND  ELECTRON  DENSITIES 

In  this  section  we  go  one  step  further  in  relating  the  intensities 
of  UV  features  to  conditions  prevailing  in  the  ionospheric  under  bombardment 
by  auroral  electrons.  We  shall  introduce  the  electron  density  n^  and  show  its 
altitude  dependence  for  the  various  incident  electron  spectra  already  discussed. 
Trere  are  certain  caveats  we  wish  to  address  before  this  is  done.  Many  para¬ 
meters  effect  ng  through  chemistry,  transport  3nd  auroral  turbulence.  Tc  avoio 
the  complicating  effect  of  plasma  transport,  we  limit  our  considerations  to  the 
E-region.  We  also  consider  situations  where  turbulence  is  weak  enough  not  to 
measurably  affect  the  distribution  of  the  plasma. 

A  range  of  calculated  ng  profiles  is  possible  with  regard  to  the 
chemistry  through  variations  in  densities  with  long  resident  times  and  un¬ 
certainties  in  rate  coefficients  which  include  those  terms  affecting  tempera¬ 
tures  as  well  as  chemically  derived  densities.  Examples  of  species  with  long 
resident  times  whose  variations  from  aurora  to  aurora  must  be  kept  in  mind 
are  NO  and  0^.  Variations  in  their  densities  can  significantly  affect  ion 
concentrations  although  produce  much  less  affect  on  their  sum  or  equivalently,  nf. 
The  concentration  of  NO  will  slowly  build  up  under  continuous  auroral  activity 
and  thus  can  depend  significantly  on  conditions  minutes  to  hours  earlier  than 
the  given  time  of  interest.  Table  1  gives  the  NO  density  as  well  as  those  .f 
the  dominant  neutral  spacies  used  to  generate  the  n£  profiles  to  follow.  The 
NO  density  comes  from  a  previous  calculation  for  modest  energy  deposition  con¬ 
tinuing  for  about  I  hour.  Although  wn  have  used  this  NO  density  profile  to 
calculate  the  ng  densities  to  follow,  use  of  substantially  different  profiles 
(factor  of  ten  either  way,  for  example)  would  not  noticeably  affect  n  in 
comparison  to  the  degree  of  effect  caused  by  changes  in  the  electron  spectra 
over  the  chosen  range  ot  spectra,  hardness. 

We  have  calculated  electron  densities  usino  a  detailed 
chemistry  model  as  well  as  the  approximation 

"e(z)  =  [-i (z)/o(r)]1/2  cm-3 
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where  ■■  is  the  total  ionization  rate  and  n  is  the  effective  recombination 
coefficient.  Densities  to  follow  come  from  the  detailed  calculations  but  do 
not  differ  significantly  from  those  obtained  using  Equation  (11).  The  re¬ 
spective  profiles  are  essentially  the  same  below  their  maxima  with  n-  30/  de¬ 
parture  at  higher  altitudes.  The  larger  values  in  the  departure  region  refer 
to  the  simplier  calculation. 

In  our  chemistry  model,  we  solve  a  set  of  coupled  rate 
equations  with  each  equation  given  by 


dn-(z)  ,  , 

— s  Pi (z) -  1  i  (z )  [cm  -s  ]  (12) 

where  p.  and  1^  are  the  production  and  loss  terms  for  the  ion  chemical  species. 

We  currently  model  16  species  which  include  metastable  as  well  as  ground  state 
ions  and  minor  neutral  species.  The  important  ion  species  are  N0+,  Cu+,  0*,  and 
N^+  with  the  latter  providing  an  important  source  to  the  others  through  charge 
transfer  reactions.  The  electron  density  is  specified  by  requiring  charge  neutrality. 
Most  of  our  rate  coefficients  come  from  Roble  and  Rees  (1977)  with  some  adjustments 
based  on  Torr  and  Terr  (1982).  A  number  of  these  coefficients  are  allowed  to  vary 
with  T  ,  the  e^ctror  temperature.  For  results  to  follow,  Tg  has  been  specified 
fr0rr  those  calculations  by  Roble  and  Rees  (1977)  which  apply  in  an  appropriate  wa\ 
to  our  conditions. 

Having  briefly  described  the  conditions  under  which  n  is  specified, 
we  now  turn  to  our  final  results  which  appear  in  Figures  13  and  14  Shown  are 
total  ionization  rates  and  corresponding  ng  profiles  for  Maxwellian  and  Gaussian 
distributions.  As  before,  energy  content  is  fixed  at  1  erg/cm2-s.  The  chemistry 
calculations  were  extended  to  sufficiently  long  times  that  steady  state  conditions 
prevail.  Our  interest  here  is  not  so  much  in  the  detailed  behavior  of  the  n  profiles 
shown  but  rather  in  their  relative  sensitivity  to  En  over  a  range  which  gives 
measurable  variations  in  the  intensities  under  consideration.  We  will,  however,  max' 
some  comparisons  between  no  profiles  in  further  discussion  below  for  the  two  chosen 
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we  do  re  >  c  :  e 


issue  of  just  how  accurately  profiles  can  be  deter- 


source  representations.  Giver  ou*  emphasis.  we  do  rot  •  c  i  se  *  t he 

5  ever  r,rtsert 

uncertainties  in  rate  coefficients  and  the  startina  condt'ors  of  tne 
ionosphere.  That  issue  must  eventually  be  addressed  to  help  answer  the 
question  of  the  practicabi 1 ity  of  using  satellite  observed  optical  emissions 
for  ionospheric  monitoring.  Returning  to  Figures  13  and  14,  it  is  dear 
that  np  cannot  be  very  well  specified  from  our  chosen  set  of  UV  features 
for  soft  electron  precipitation  such  as  given  by  Em  in  the  range  .5  to  say 
1.5  keV.  Over  this  range,  intensity  ratios  change  modestly  whereas  significant, 
variations  occur  in  the  distribution  of  r,  For  harder  electron  spectra,  a 
more  precise  specification  of  np  is  possible  due  both  to  the  stronger  intensity 
variations  as  shown  in  Figures  7  and  8  and  to  the  more  confined  nature  of 
the  ng  profiles  about  their  maxima. 


In  the  previous  section,  we  discussed  tre  pcss:  t  ,  1 1 t.<  o*  direct'., 
estimating  the  altitude  of  maximum  energy  deposition  zn  from  the  chosen  inten- 
sities.  Figure  10  showed  that  such  a  specification  of  zp  improves  as  it  de¬ 
creases  in  value  and  that  below  a  value  of  •u  120  km,  the  techmoue  appears 
attractive.  For  such  values,  zp  also  applies  to  ng  which  thus  allows  us  tc 
directly  estimate  the  altitude  of  maximum  n0  from  results  like  those  ir  Cicure 
10.  The  fact  that  the  energy  deposition  rate  and  np  peak  at  the  same  altitude 
(provided  that  the  electron  spectrum  is  not  too  soft)  can  be  seen  by  comparing 
the  panels  in  either  Figure  13  or  14.  It  should  be  added  that  the  ion  produc¬ 
tion  rates  shown  have  the  same  shape  as  the  energy  deposition  rate  to  which 
Figure  10  refers. 


To  conclude  this  section,  we  briefly  comment  on  c^e'-erce s 

between  n  profiles  for  the  two  source  representations  in  which  respective 
e 

characteristic  energies  are  chosen  to  give  the  same  altitude  for  maximum 

ionization.  We  have  in  mind  the  situation  where  the  altitude  of  np  is 

,  ma  x 

approximately  known  using  information  such  as  given  in  Figure  iO.  i he 

question  then  arises  as  to  how  much  variation  can  be  expected  in  nf  over 
the  possible  range  of  incident  electron  spectra  which  give  the  same  altitude 
for  n^  .  In  figures  13  and  14,  we  see  that  E^  =  2.5  keV  and  Eg  =  7.5  keV 
lead  toaa  common  altitude  of  n  10?  km  for  ng  and  thus,  let  us  address  these 
energies.  As  expected,  ng  is  more  broadly  dTllr ibuted  for  the  Maxwellian 
source  and  accordingly  has  a  smaller  maximum.  The  difference  in  effective 
lower  boundaries,  as  defined  here  to  be  where  n?  is  falling  rapidly,  is  4  km. 


1 


I 


4 


in 


At  higher  altitudes,  such  as  160  km,  ng  for  the  Maxwellian  source  is  a  facto 
of  v  1.5  larger  in  magnitude.  Such  differences  appear  modest  under  the  con¬ 
ditions  for  which  they  arise.  Nevertheless,  whether  such  differences  are 
important  depends,  of  course,  on  how  one  wishes  to  use  such  deduced  profiles 
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Tne>'e  are  rany  mechanisms  one  must  consider  in  relating  aurora' 
optica’  properties  to  associated  plasma  conditions.  Those  considered  in  this 
paper  have  been  1}  electron  transport  based  on  a  particle-particle  description, 
2)  excitation  and  ionization  associated  with  this  transport,  3)  pure  absorption, 
ana  4)  the  chemistry  which  relates  excitation  and  ionization  rates  to  the  elec¬ 
tron  and  ion  densities.  Other  mechanisms  identified  were  photon  scattering  ana 
plasma  turbulence.  Scattering  was  noted  to  have  no  effect  on  the  features  con¬ 
sidered  here  with  the  exception  of  01  1356  A  for  which  minor  effects  occur. 

This  statement  is  dependent  on  viewing  conditions  and  here  we  are  referring  to 
nor-hori zontal  loot  directions.  Turbulence  may  manifest  itself  in  modifying 
energy  deposition  of  the  incoming  energetic  auroral  electrons  or  in  creating 
i rregul ar i t ies  in  the  plasma.  Since  neither  effect  has  been  addressed  here, 
our  conclusions  relate  to  optical  properties  based  on  purely  particle-particle 
interactions  and  a  stable  plasma. 

Within,  the  above  confines,  there  are  various  issues  to  be  addressed 
which  we  have  begun  to  do  in  this  work.  They  relate  to  1)  dependence  of 
emissions  on  the  neutral  atmospheric  composition  and  2)  sensitivity  of  these 
emissions  to  the  incident  electron  spectrum.  One  would  hope  that  a  measure  of 
the  richt  set  of  intensities  could  yield  both  composition  information  and  a 
useful  description  of  the  incident  electron  spectrum.  By  useful,  we  me  a-  tnat 
the  descriotior  is  adequate  to  allow  specification  of,  say.  the  associated 
electron  density  altitude  profile  to  within  some  limit  of  tolerance  decree  p, 
its  appl ication. 

We  have  not  yet  carried  out  our  analysis  to  the  point  of  providing 
firm  answers  to  questions  of  a  practical  nature  related  to  degrees  of  uncerta ' nty 
in  the  deduced  parameters.  It  is  clear,  though,  that  accurate  measurements  of 
intensities  such  as  those  considered  here  will  have  useful  diagnostic  value  m 
establishing  dominant  regions  of  excitation  and  ionization. 
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5.6  (13) 

1.5  (13) 
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3.2  (14) 

1.0  (14) 

3.2  (10) 

'aMe  1.  Model  etmosDhere  used  in  the  calculations.  The 

N^,  C  ^ ,  and  0  densities  come  from  Jacchia  il9"’7). 
The  NO  density  comes  from  chemistry  mode-lire  by 
us. 
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"1 

v ' ,  v" 

EMISSION 

02  PURE  ABSORPTION  1 

(A) 

FACTOR  (%) 

CROSS  SECTION  (10*16cr/; 

1273 

6,0 

1.5 

G .  09 

1325 

4,0 

4.9 

2.0 

1354 

3,0 

6.0 

7.1 

(l  384  | 

(  2,0  j 

(4.9  | 

13.0 

jl  382  j 

I5’2! 

V-°\ 

3914 

0,0 

64.8 

- J 

0.0 

Table  3  Emission  factors  and  0^  pore  absorption  cross  sections 
for  the  3914  A  and  LBH  bands  under  investigation. 

The  emission  factors  are  given  as  percents  of  the  total 
a  production. 
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FIGURE  CA 


Figure  1. 


Figure  2. 


Figure  3. 


Figure  4. 


Figure  5. 


Figure  6 . 
Figure  7. 


Figure  8. 


Figure  9. 


*  Jil  , 


An  Auroral  FUV  spectrurr  from  P.D.  Feldman  (private  communica¬ 
tion,  1981)  obtained  in  a  rocket  experiment  with  the  spectrometer 
looking  down  from  above  the  emitting  region.  Dashed  curve  gives 
the  photo  absorption  cross  section.  Altitude  scale  on  right 
is  to  be  used  with  adjacent  cross  section  scale. 

Nadir  viewing  FUV  photometer  data  obtained  on  satellite  S  3-4 
(Huffman  et  al.,  1981).  Peak  sensitivity  occurs  at  1550  A. 

Data  show  latitudinal  distributions  of  the  continuous  and 
discrete  auroras. 

0^  pure  absorption  optical  depths  for  the  LBH  bancs  under  ir, edi¬ 
tion.  These  are  obtained  using  the  Jacchia  (1977)  model  atncsohe 

Density  ratio  profiles  with  the  ratio  given  by  the  density 
values  from  Table  3  (Jacchia,  1977)  divided  by  the  corresponci^c 
values  from  the  model  reported  by  Rees  et  al .  (1977  )  and  She’-: 
et  al.  (1979). 

N2  ( 4 ^7! g )  volume  production  rates  for  incident  electron  spect-a 

given  by  Maxwellian  (upper  panel)  and  Gaussian  (lower  pane1) 

2 

distributions.  All  incident  fluxes  contain  1  erg/cm  -s. 

Similar  to  Figure  5  except  for  0(  S). 

Nadir  viewing  intensities  versus  characteri stic  energy  E^  for 

incident  electron  spectra  given  by  Maxwel’ian  di stributions . 

2 

Each  distribution  contains  1  erg/cm  -s. 

Similar  to  Figure  7  except  for  incident  electron  spectra  given 
by  Gaussian  distributions. 

Nadir  viewing  intensities  versus  incident  electron  energy 
content  Qq  for  Maxwellian  distributions  with  E^  =  2.5  keV. 


Figure  10.  Nadir  viewing  intensities  versus  Zp,  the  altitude  of  peak 
energy  deposition. 

Figure  11.  01  1356  A  and  LBH  1354  A  nadir  viewing  intensities  for  two 

n ( 0 )  models. 

Figure  12.  LBH  nadir  viewing  intensities  for  three  different  C>2  density 
distributions.  The  unsealed  distribution  is  that  given  in 
Table  1. 

Figure  13.  Total  ion  production  rates  and  electron  densities  for  Maxwellian 
distributions  with  =  1  erg/cm-s. 

Figure  14.  Similar  to  Figure  13  except  for  Gaussian  distributions. 
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